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We show that a coupling between chameleon-like scalar fields and photons induces linear and 
circular polarization in the light from astrophysical sources. In this context chameleon-like scalar 
fields includes those of the Olive-Pospelov (OP) model, which describes a varying fine structure 
constant. We determine the form of this polarization numerically and give analytic expressions 
in two useful limits. By comparing the predicted signal with current observations we are able to 
improve the constraints on the chameleon-photon coupling and the coupling in the OP model by 
over two orders of magnitude. It is argued that, if observed, the distinctive form of the chameleon 
induced circular polarization would represent a smoking gun for the presence of a chameleon. We 
also report a tentative statistical detection of a chameleon-like scalar field from observations of 
starlight polarization in our galaxy. 
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I. INTRODUCTION 

Extensions of the Standard Model of particle physics, 
such as string theory, introduce many new scalar fields 
which are not seen in the Standard Model. Such scalar 
fields are commonly invoked to explain the observed ac- 
celeration of the universe, as inflation [l| or dark energy 
2| fields, or to cause variations in fundamental constants 
^. If new scalar fields do indeed exist in the Universe, it 
is important to understand the properties of the theoret- 
ical models that describe them, e.g. the interactions of 
the scalar fields with themselves and with matter, which 
may give rise to additional observable effects that could 
be tested and constrained by experiments. 

In this article we consider the effect of scalar field theo- 
ries with a self-interaction potential, V{(j)), and couplings 
to matter and light on observations of the polarization of 
light from astrophysical sources. These scalar field theo- 
ries are described by the action: 




%^-V{^) (1) 



where the Sm are the matter actions for the matter fields 
1pm , and the functions Bi{(f>/MQ) and Bf{4>/AIo) deter- 
mine the couplings of the scalar field, 4>, to the i^^ matter 
species, ipi, and to the photon field respectively. A scalar 
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field with such couplings to matter fields might be ex- 
pected to give rise to fifth force effects or violations of 
the weak equivalence principle. In this article we are 
specifically interested in a scalar field, (j), which is light 
in relatively low density regions such as galaxies, galaxy 
clusters and the inter-galactic medium. More precisely, 
in these regions we require the mass, m<^, of small per- 
turbations about the background value of the scalar field, 
(j)h, to satisfy < lO^^^eV/c^. Hence, any force me- 
diated by 4> would have a range = l/m^ > 20 km. 
Additionally we require that coupling between photons 
and (p in these regions is relatively strong: 



> 10-" GeV 
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Therefore even if the coupling to matter is much weaker 
than the coupling to photons, since roughly 10~* of the 
mass of nuclei is due to electromagnetic interactions the 
^-mediated force between individual nuclei in these back- 
grounds will be at least 10^ times the strength of gravity 
on scales smaller than A^. 

One might, at first glance, conclude that a scalar field 
theory with these properties is already strongly ruled out 
by laboratory constraints, e.g. [1, [E S S H, HI, on the 
strength of fifth forces. Specifically, measurements of the 
displacement of a micro-machined silicon cantilever us- 
ing a fibre interferometer, reported in require that 
a Yukawa-type fifth force with strength 10'' times that 
of gravity has a range A^ < 5 fim with 95% confidence. 
This, however, does not rule out the models we wish to 
consider as neither the strength nor the range of the 0- 
mediated force are necessarily the same in the relatively 
high density environment of the laboratory as they are 
in the low density background of space. 

In recent years, two classes of models have arisen that 
allow a scalar field that is strongly interacting in low 



2 



density environments and yet is currently undetected in 
laboratory tests: the chameleon model, [13) [HI) and the 
Olive-Pospelov (OP) model [H. Both models are de- 
scribed in detail in the following Section. The mecha- 
nism by which these models avoid laboratory tests can 
be understood by extremizing Eq. ([1]) with respect to (j) 
to give the following field equation: 



□0 = K=ff,0(</')T,,,FV4) 



where 



Kff(0;T„F2) = T/(0) 



(2) 



Tm, (3) 



and Tin = 9iJ.vT^-^ is the trace of the energy momentum 
tensor for matter: 



2 5S^ 



-Pm, where pm is 
Both the chameleon 



For non-relativistic matter : 
the energy density of the matter 
and OP models play the scalar field potential, off 
against the matter couplings. Bp and B^, to make the 
vacuum expectation value (VEV), and hence the prop- 
erties of the field, depend strongly on the local den- 
sity of matter. For convenience we shall refer to such 
a scalar field (f) as the chameleon or chameleon field, 
however our analysis applies equally well to both the 
chameleon and OP models. In this analysis we posit 
a universal coupling to the different matter species i.e. 
Bi{(j)/Mo) — Bm{(j)/AIo). Although it is not required by 
either model, we make this assumption because it simpli- 
fies the analysis whilst having little effect our conclusions. 
The best constraints on Mq come from the requirement 
that corrections to particle physics are small, which lim- 
its Mo > lO^GeV [i3. 

A coupling between matter fields and the chameleon 
potentially causes violations of the weak equivalence 
principle (WEP) and other fifth force effects such as an 
effective alteration to Newton's inverse square law. A 
coupling between photons and chameleons introduces ad- 
ditional observable phenomena for the chameleon field. 
If such a coupling has super-gravitational strength, it 
can result in a non-negligible conversion of photons to 
chameleons and vice versa. The detectable effects asso- 
ciated with this conversion are similar to those predicted 
for axion-like-particles (ALPs) which interact with light 
[T^ . Mixing requires the interaction between two pho- 
tons and one scalar particle, and so the effects of the 
mixing are most likely to be seen when a photon, or 
a scalar particle is passing through an external electro- 
magnetic field. The chameleon-photon coupling induces 
both birefringence and dichroism [Tsl . [T4 | in a coherent 
photon beam passing through an external magnetic field. 
These effects could be detected by laboratory searches, 
such as the polarization experiments PVLAS, Q&A, and 
BMV [ll,[li,[i7|,|]J|, that are sensitive to new hypothet- 
ical particles with a small mass and coupling to photons. 



Such experiments can constrain the coupling M in the 
chameleon model, indeed the otherwise anomalous de- 
tection of birefringence with a 5.5 T magnetic field by 
PVLAS [l^, could, at least in principle, be explained 
by the presence of a chameleon field [ij]. For the most 
widely studied class of potentials, the PVLAS data was 
found to rule out M < 2x 10*^ GeV [H. In the OP model 
the mass of the scalar field in the laboratory is too large 
to produce a detectable effect in these experiments. 

If chameleons exist and couple to photons, then they 
could, as suggested in [l^, [20|, be trapped and slowly 
converted back into photons resulting in a long lived 
chameleonic afterglow. A number of experiments, most 



notably GammeV 21| , are searching or aiming to search 
for this afterglow effect. The GammeV chameleon search 
recently announced its first results which, for models 
with < 10~'^eV in the interior of the experiment, 
ruled out 2.4 x 10^ GeV < M < 3.9 x 10^ GeV [2^. Ul- 
timately GammeV may be sensitive to M < 10^ GeV 
[Tot , and an optimal sensitivity for afterglow searches of 
M < 10^'^ GeV is feasible within the constraints of cur- 
rently available technology [19]. Indeed for any of the ef- 
fects associated with the coupling to photons to be large 
enough to be detected in the laboratory, either now or 
in the foreseeable future, one must have M < 10^" GeV. 
Such laboratory constraints do not, however, apply to 
the OP models. 

These laboratory experiments need to be performed in 
a very good approximation to a vacuum otherwise the 
chameleon becomes too heavy to have a noticeable ef- 
fect. A complimentary approach to testing chameleon- 
photon couplings is to look for the effects of the coupling 
in observations of astronomical objects. The densities of 
interstellar space are typically very low and so the effects 
of the chameleon may be significant. Light from all as- 
tronomical objects travels a significant distance through 
magnetic fields in galaxies, galaxy clusters and possibly 
in the intergalactic medium, before reaching the earth. 
However astronomical magnetic fields are typically made 
up of large numbers of randomly oriented magnetic do- 
mains - a very different scenario to the well controlled 
constant magnetic fields of laboratory experiments. In 
contrast to laboratory tests, such astrophysical effects 
should be see in the OP as well as the chameleon models. 

The coupling between photons and chameleons means 
that photon number is not conserved, however, as the flux 
of photons emitted by astronomical objects is diflicult to 
determine, measurements of flux cannot be used to bound 
the parameters of the chameleon model. In the following 
sections we show how the coupling between photons and 
chameleons generates polarization in the light from astro- 
nomical objects. Therefore measurements of polarization 
can be used to constrain the parameters of the chameleon 
model because the intrinsic polarization of astronomical 
objects is often very well constrained. Astronomers are 
interested in measuring polarization because it can pro- 
vide information both about the source of the radiation 
and about any magnetic fields present between the source 
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and the earth. Very precise astronomical polarization 
measurements are therefore available, and can be used 
to constrain the chameleon model. 

This article is organized as follows: in fJII]we introduce 
and provide further details of the two classes of model 
to which our analysis applies i.e the chameleon and OP 
models. The coupling of a chameleon-like-particle to pho- 
tons, is essentially the same as that which is assumed for 
a scalar axion-like particle (ALP). There is a great deal 
of literature concerning constraints (both local and as- 
trophysical) on ALPs. However, the density dependent 
mass of a chameleon field, allows chameleon theories to 
evade the tightest of these constraints. In mill we review 
previous constraints on ALPs, and consider to what ex- 
tent they do, or do not, apply to chameleon-like models. 
In mVl we consider how the existence of a chameleon-like 
field alters the polarization of light from astrophysical ob- 
jects as it passes through an astrophysical magnetic field, 
and derive the form of the induced polarization. In ifVl we 
discuss the observed and predicted properties of the dif- 
ferent types of large scale astrophysical magnetic fields. 
In SjVTl we apply the results of the previous sections, and 
use astrophysical polarization observations to constrain 
the chameleon to photon coupling. We find that such 
measurements place the tightest constraints yet on this 
coupling. Applying the analysis to starlight polarisation 
in our galaxy we find a tentative statistical detection of 
a chameleon-like scalar field. Finally we summarize our 
results in ^VlIIl 



II. THE MODELS 

A. Chameleon Model 

In the chameleon model the coupling functions and 
_Bf in Eq. ([1]) are well approximated by linear functions 
of (p i.e. Bp ^ 1 + 4>/M and « 1 + 20/Mo; for rea- 
sons of naturalness, if 1/M ^ 0, M ~ 0{Mo) is usually 
assumed. The strength of the matter coupling is deter- 
mined by Bp ^p and Brn,.4>, hence in the chameleon model 
the coupling strength does not depend explicitly on the 
VEV of 4>. The model was originally proposed by Khoury 
and Weltman with Mq ~ O(Mpi), which results in a 
gravitational strength coupling between matter and the 
chameleon field, </). The ability of the chameleon with this 
coupling to behave as dark energy was discussed in . 
The coupling to photons, Bp, was not expressly consid- 
ered in [10] although with M ~ Mq ^ Mpi the most 
pronounced new effect is a virtually undetectable density 
dependence in the fine structure constant. With a grav- 
itational strength coupling to matter (and possibly also 
photons) the chameleon field could be detected by labo- 
ratory, satellite, solar system and astrophysical tests (e.g. 
structure formation [2J]) of gravity. A potentially much 
wider phenomenology was opened up, when Mota and 
Shaw [ll| showed that coupling between chameleon fields 
and matter could be many orders of magnitude stronger 



than gravity, Mq <C Mpi, and yet still be compatible with 
all existing experimental data. The properties of such 
strongly coupled chameleon fields can probed u sing ex- 
periments designed to measure the Casimir force jlll. [25| . 
In addition, with a strong coupling and M ^ 0{Mo), 
Brax et al. noted that interactions between 

chameleons fields and photons would result in potentially 
detectable effects similar to those predicted for axion- 
like-particles (ALPs) which interact with light. 

It should be noted, that it is generally seen as 'natural', 
from the point of view of string theory, to have M « Mpi. 
This relation also arises in /(i?) modified gravity theo- 
ries (see e.g. Ref. [l^ and references therein). It has 
also been suggested, however, that the chameleon field 
arises from the compactification of extra dimensions, [2^ . 
In this case, there is no particular reason why the true 
Planck scale (i.e. that of the whole of space time in- 
cluding the extra-dimensions) should be the same as the 
effective 4-dimensional Planck scale defined by Mpj. In- 
deed having the true Planck scale much lower than Mpi 
has been suggested as a means of solving the Hierar- 
chy problem (e.g. the ADD scenario [1^). In string- 
theory too, there is no particular reason why the string- 
scale should be the same as the effective four-dimensional 
Planck scale. It is also possible that the chameleon might 
arise as a result of new physics with an associated en- 
ergy scale greater than the electroweak scale but much 
less than A'/pi. Therefore in this article we consider M 
as a free energy scale to be constrained by experiment. 
This said, to date, no one has managed to find such a 
chameleon theory (with either M ~ Mpi or otherwise) in 
the low-energy limit of a more fundamental high energy 
theory (e.g. supergravity). 

The chameleon model evades the strong constraints 
imposed by local tests of gravity [13, [HI through non- 
linear self interactions of the field described by the po- 
tential T^(0), and hence the field may couple with super- 
gravitational strength. 'Non-linear' in this case means 
that V,fj, is a strongly non-linear function of 0, and the 
mass of the scalar field, = \/V^^^{4>), therefore de- 
pends strongly on the VEV of 0. The VEV of in a 
given background is determined by the minimum of the 
effective potential, and therefore the position of the 
minimum depends on and . V{4>) is chosen so that 
is larger in high density regions than it is in low den- 
sity regions. It is then possible to ensure that in a galaxy 
> 20 km whereas in the laboratory < 5 /im. As- 
suming that InBp ^j, > and ln_B„j_0 > 0, the chameleon 
mechanism requires: 

< 0, V;^^ > 0, ^000 < 0. 

To provide some intuition for what we expect to be 
in a low density region such as a galaxy or galaxy cluster, 
we consider the most widely studied class of chameleon 
models where: 



V{4i) « const. -I- 
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with 0/A < 1, n > -1 and n 0{1). We note that 
this includes potentials with the form V — const. — 
A"* ln((/)/Af). A is constrained by experiments to be at 
most a few orders of magnitude larger than the dark 

energy energy scale Aq = p^^"* = (2.4 ±0.3) x 10"^ eV 
[ill [23] . When the chameleon is posited as an explana- 
tion for dark energy, it is therefore considered natural 
to take A « Ao [10|, lill- The minimum of the effective 
potential occurs when (j) — (j>i, where — V^0(0(,) — Pb/M^. 
The mass of the chameleon at this minimum is given by 



VK^^I^fJ, so: 



— A\/n + 



Pb 



n + 2 
2(" + l) 



Assuming Mq « M > 10^ GeV (i.e. the region which 
is not currently accessible to laboratory experiments), 
we have < 10~^^eV for all —l<n< 5.6 in a 
background, such as a galaxy or galaxy cluster, with 

Pb « 10- 



^''gcm 



B. Olive-Pospelov Model 

The Olive-Pospelov (OP) model [l^ was proposed as 
a way to allow particle masses and coupling 'constants' 
to depend on the local energy density of matter. The 
model could therefore provide an explanation for the 
6(7 difference between the value of the fine structure 
constant, a = e^/hc, in the laboratory and that ex- 
trapolated from the spectra of 128 QSO absorption sys - 
tems at redshifts 0.5 < 2: < 3 by Webb et al. [29j : 
Aa/a = (aqso-aiab)/aiab = -0.57±0.10x 10"^. There 
is now a great deal of tension between other potential the- 
oretical explanations for this data see e.g. Refs. (sl. [soj. 
and the most recent local atomic clock constraints on any 
local time variation of a [sil [3^ . Issj . In the OP model, a 
is locally time-independent and hence these constraints 
are avoided. 

The OP model could also describe a density depen- 
dent electron-proton mass ratio fi = nip/ me- Rein- 
hold et al. ,34j reported a 4(T indication of a vari- 
ation in /i. They analysed the H2 wavelengths of the 
spectra of two absorbers at z « 2.6 and z « 3.0 ob- 
served using the Very Large Telescope (VLT), finding 
Afi/jj. — 24.4 ± 5.9 X 10~^. It was shown, however, in 
Ref. |i35i] that due to wavelength calibration errors in the 
spectrograph on the VLT identified in Ref. [s^, the re- 
sult of Reinhold et al. could no longer be trusted. The 
reanalysis performed by King et al. [ssj , in which data 
from an additional object at z w 2.8 was also included, 
found Afi/fi = (2.6 ± 3.0) x 10~®, which is consistent 
with no change. Very recently Levshakov et al. (37j 
have reported evidence for a spatial variation in /x found 
by measuring ammonia emission lines in the Milky Way: 
Sp/p = -(4-14) X 10-^ 

All of these astronomical measurements of fj, and a, 
were made in regions where the average density of matter. 



Pb is very low compared to the ambient density of matter 
in a laboratory. The background density for all of these 
measurements, pb, is similar to the average density of a 
galaxy or galaxy cluster i.e. pb ^ 10~^^gcm~'^. These 
measurements could therefore be an indication that some 
or all of the 'constants' of Nature depend on the ambient 
density of matter. The OP model realises just such a 
density dependent variation in a manner that does not 
conflict with local tests of gravity. In this model, the 
coupling functions, Bm and Bp, are chosen so that they 
are close to their minimum (which occurs at (j^m)- 
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Bf = 1 + ^ 



Mo 



1 



Mo 



For reasons of naturalness, one would expect ^F,6n ^ 
0(1) [l2|- In contrast to the chameleon model, the OP 
model does not require that the potential, V{<j)), con- 
tain non- linear self interaction terms, and in the simplest 
model: 



Vi<t>) 



At 



In a background with density pb, assuming <C Pb, 
as is usually the case, and fixing the definition of Mq by 
setting = 1, the value of (j) at the minimum of the 
effective potential, ^min is given by: 



Ph 



Pb + Af 



(4) 



In the laboratory environment, pb 3> A|, and so (ftmin ~ 
(/),„. Additionally the effective matter coupling is small 
enough to evade experimental constraints. In low density 
regions such as galaxies Af ^ pb, so (pmin ~ 0. The 
change in a between the laboratory and a low density 
region such as galaxy is given by: 



Sa 
a 



aiab 



"lab 



2 



Mo 



To explain the Webb et al. value [29i] of Aa/a, one 

would require 4>m/Mo ~_3£p x 10~^. 

Olive and Pospelov [IJl found that the current best 

constraints on Mo are Mq > 15TeV and Mo^^^^^ > 
3TeV. We define rrC'^'' = Aj/Mo to be the mass of small 
perturbations in (j> in a low density region (i.e. ph <C Af), 
and let A™'^ — 1 /m^^ specify the range of the (/)-mediatcd 
force in such a region. It was found in [l^ that 



10-3M, 



M, 







ITcV 



1 km 



< 10-^ 



10" 



For there to be measurable differences between the parti- 
cle masses and coupling constants in the laboratory and 
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in regions with p « 10 ^^gcm ^, one must require: 



Mo \ lkm\ ^ 7 



In the low density regions where w the effective 
coupling to the photon field for small perturbations in 
is: 



1/M : 



dln^F 




-10"^ GeV 



It is clear then that a field with the required properties, 
5077 ^ 10~" GeV^^ and A™'' > 2 km, is perfectly com- 
patible with current experimental constraints. 

We note that in the value suggested for Ai, 

which is compatible with all current constraints, is Ai ~ 
0(l)eV. Now: 



TO0 



VleVy V Mo ) ' 



so for Mo > 15 TeV, we have < 7 x 10 eV which 
corresponds roughly to > 2800 km. 



III. CONSTRAINTS ON 
AXION-LIKE-PARTICLES 

Axion-like-particles (ALPs) can either be scalar or 
pseudo-scalar fields which couple to the electromagnetic 
field strength. If it were not for the chameleon mech- 
anism (i.e. the density-dependent mass) present in the 
chameleon and OP models, they would essentially de- 
scribe a standard scalar ALP. There are tight constraints 
on the coupling, g^^^, of ALPs to photons. In the pre- 
vious section, we discussed constraints from local experi- 
ments on chameleon-like particles, however such particles 
are also constrained by searches for ALPs. For a recent 
review of the astrophysical constraints on ALPs see Ref. 
[3^ and reference therein. In all cases, these constraints 
only apply when the ALP mass, m^, lies within a certain 
range e.g miow < < mhigh- The mass of a chameleon- 
like particle is not, however, fixed and so applying these 
constraints to chameleon models is non-trivial. We must 
take great care to identify the ambient density of the re- 
gion wherein the constraint on is required to derive 
the bound on g^j-y. 

The strongest astrophysical constraints in Ref. [38| 
come from axion production in the cores of stars. The 
application of the constraints of solar axion production 
to chameleon-like models has previously been studied in 
Ref. [l^l and [3§]. The Sun may be a powerful source 
of ALP flux, and the predicted effects of the loss of en- 
ergy of the Sun through ALP emission allows one to con- 
strain the coupling g^-yy. It must be noted, that all solar 



ALP constraints require that the ALPs actually escape 
the sun. The strongest solar ALP constraints come from 
limits on the solar neutrino flux, and this gives: 

5077 ^ 5 X 10~^"GeV"\ 

Similar constraints result from the CERN Axion Solar 
Telescope (CAST) which attempts to directly detect so- 
lar axions. However it was shown in Ref. [iJl, that solar 
chameleon-like ALPs would generally bounce off, rather 
than enter the CAST instrument, and so the CAST con- 
straints cannot be applied to chameleon models. Similar 
constraints are found from the life-time of Helium burn- 
ing (HB) stars in globular clusters: gip^^ < 10-i°GeV"^ 
Solar axion constraints are derived from production of 
axions in the solar core by the Primakoff process. In 
this region the temperature is T « 1.3 keV, and the typ- 
ical density is 150 gcm^'^. In the Helium burning stars, 
T sa lOkeV and p w 10'*gcm~'^. All of the constraints 
assume that ^ T. It was shown in Ref. [s^, that all 
solar axion production bounds are evaded if > lOkeV 
in the solar core. Similarly, the Helium burning star con- 
straints are effectively evaded if > 30keV in their 
cores. For example, if one considers a chameleon po- 
tential like k'^{k/4>), where A « 2.3 x lO-^eV, one 
finds that with a matter coupling of 10^" GeV, we have 
TO^ w 1.5MeV in the core regions of the HB stars, and 
K, 64 keV in the solar core. With the different choice 
of potential, A^ exp(A/(/)), one finds w 30keV in the 
cores of Helium burning stars but « 3 keV in the Sun 
when M = 10^° GeV. Thus the solar and HB star axion 
constraints on g^-yy = 1 /M would apply to the latter po- 
tential with M « 10^° GeV but be evaded by the former. 
If we took M « 2 x 10^ GeV, then both potentials would 
evade these constraints. Thus in the chameleon model, 
whether or not these astrophysical constraints are rele- 
vant depends greatly on the properties of the potential, 
and in particularly how it determines the behaviour of 
the theory at high densities. In general, these properties 
cannot be inferred from the low-density behaviour of the 
theories. We are concerned only with the low density 
behaviour in this work. 

At high densities, the scalar field in the OP model 
couples quadratically, rather than linearly (as an ALP 
would) to the QED term. In this way it avoids astro- 
physical constraints related to axion production in high 
density regions. 

Recently, in Ref. [i^, it was shown that polarization 
measurements of 7-ray bursts could be used to constrain 
axion production at the source of the burst. Whilst later 
in this article we will consider the potential constraints on 
chameleon-like fields from 7-ray burst polarization mea- 
surements, we will be interested in constraining any po- 
larization that is induced by the chameleon as the light 
from the 7-ray burst passes through low-density mag- 
netized regions of space (e.g. the interstellar and inter- 
galactic mediums). We will assume that axion produc- 
tion in the immediate vicinity of the 7-ray burst itself is 
negligible. In Ref. [40], it is assumed that, in the vicin- 
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ity of a 7-ray burst, there is a magnetic field of strength 
B ~ 10^ G over a distance of about ^grb ^ 10^ cm. 
Given this, it is found that: 

5077 ^ 5 X 10"^2GeV"\ 

for 8 X 10~^ eV < < 3.5 x 10^^ eV. For larger values 

of 7710 : 

g^,, < 2.2 X 10-« (^) GeV-i. 

It is noted that in the vicinity of the GRB, 77e w 
10^*^ cm~'^, corresponding to sa 2 x 10~^'*g. The effec- 
tive 'energy density' to which the chameleon field couples 
is not just but ptot = Pm + — E'^/2. Thus for 

the GRB, ptot « « 4.4 x lO^^gcm-^. Such a ptot 

places one in the high-density region of the OP model, 
where the only couples to photons quadratically and 
hence no longer behaves as an axion. In chameleon the- 
ories, if ¥{(/)) = AV(0/A) where /'(I) - /"(I) ~ 0(1) 
and A « 2.3 x 10~^ eV (as is usually assumed), one finds 
that 7710 ^ 10"'^ eV when 

ptot » MA' « 2.8 X 10-«gcm-3 (^qT^) . 

Thus the strongest constraint on g^^j from Ref. [i^ 
does not apply here. If we take V — A'' (A/0) or 
V — A'*exp(A/0) then we predict 771^ « 0.8 eV or 
« 0.4 eV respectively for M « lO^^GeV; hence 
M « 10^° GeV is allowed. Indeed, we find that the bound 
of Ref. [4^ , would allow all such chameleon models with 
M > 10^ GeV. 

It should also be noted, that axion-like chameleon 
production from the magnetic fields of neutron stars 
would also be greatly suppressed. For a neutron star 
B w 10^^ G, which corresponds to ptot ~ 44gcm~'^ and 
hence a very heavy chameleon particle. 

It is clear then that astrophysical ALP constraints 
coming from relatively high density regions do not ap- 
ply to the OP model, and the extent to which they apply 
to a chameleon theory depends greatly on the precise 
choice of potential. For at least one popular choice of 
potential {V = A'*(A/0)) one of the constraints noted 
above applies. Furthermore, because the chameleon field 
is very heavy in high density regions, we expect any ini- 
tial chameleon flux from stars or objects to be greatly 
suppressed relative to that which one would expect for a 
standard ALP. 

There has also been a great deal of work on conversion 
of photons to very light ALPs in relatively low density 
backgrounds ( e.g. the interstellar medium). See, for ex- 
ample, Refs. [41I ,42 . .43. , 44, 45, 46], and for a recent 
review see Ref. [43|- In relatively low density regions, 
chameleon-like particles behave essentially like standard 
axion-like-particles. Therefore much of the analysis pre- 
sented in the aforementioned works is directly applica- 
ble. Only where a initial axion flux from, for example. 



a star or quasar has been assumed will the analysis dif- 
fer. Many of these studies have focused on photon-axion 
conversion in the inter-galactic medium. Magnetic fields 
with strength B ~ 10~^ G are generally seen as plausible 
in the inter-galactic medium. It is suspected that such 
fields would be coherent over scales of about a mega- 
parsec or so. We discuss this further in S|Vl For rea- 
sonable values of the electron number density, 77c, in 
the inter-galactic medium, it is commonly found that 
5077 < 10"^"GeV(lnG/BiGM) or so 143]. 

Carlson and Garretson [i^] specifically considered the 
effects of photon to ALP conversion induced by the mag- 
netic field of our own galaxy. This discussion is directly 
relevant to our work. In their work they were only able 
to constrain g^^-y < 10~^ GeV, however they suggested 
a method that would allow couplings down to 10~^ GeV 
to be probed. In our work to use a different method 
to constrain g^^^ down to 10~^ GeV. Ref. [3] is also 
interesting because it is noted that small scale fluctua- 
tions in the electron-density can lead to an enhancement 
of the photon to ALP conversion rate. In their work, 
the enhancement effect was estimated to be very large 
for visible light. We discuss this further in Appendix [X[ 
and note that the size of the enhancement effect found in 
Ref. [3] was in part due to, what is now, an old model 
for the electron-density fluctuations. Using the more re- 
cent NE2001 model [4^, we show in Appendix [51 that 
the enhancement effect is expected to be no larger than 
0(1) in the local interstellar medium. Due to the com- 
plexities and additional uncertainties associated with the 
structure of electron-density fluctuations, particularly at 
parsec scales, which determine the magnitude of any en- 
hancement, we have neglected the potential enhancement 
effect of Ref. [13] from our analysis. As we note in Ap- 
pendix]^ however, we do not expect this to greatly alter 
our conclusions. Similarly, the analysis of Ref. [1^ is 
applicable to chameleon-like fields, however our analysis 
goes beyond what was presented there. 

We also comment on Ref. [4^. Here a supercluster 
magnetic field with strength 1 pG coherent over a scale 
of 10 Mpc was assumed. Additionally an enhancement ef- 
fect similar to that derived in Ref. [4J] was employed. It 
must be noted that the magnitude of any enhancement 
effect depends greatly on both the magnitude and the 
spatial scale of the spectrum of electron-density fluctua- 
tions. The former is fairly well known for electrons in our 
galaxy, whereas the latter is less well known. In the con- 
text of electrons in a supercluster neither is well known. 
Additionally evidence for a field strength of i? w 1 /iG co- 
herent over 10 Mpc was tentative at best at the time that 
Ref. [4^ , and a more recent analysis [i^ suggests that if 
such a field docs exist it is either weaker, B ~ 0.1 pG, or 
only coherent over much smaller scales ~ 100 kpc. Even 
if such a field does exist, it is also not clear precisely what 
distance along the line of sight the field extends. As such, 
the constraint: gcp^j < 10-13 GeV~^ quoted in Ref. [H 
relies on many assumptions, with at best only tentative 
observational support. Removing any one of these as- 
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sumptions, would allow for much larger couplings. In 
this work, we are primarily concerned with constraints 
on photon to chameleon conversion in astrophysical re- 
gions where there is strong evidence for magnetic fields, 
and the properties of such magnetic fields are relatively 
well known. We also note in Appendix El that making 

1 /2 

the reasonable assumption <^((5rie)^) / ("-e) ~ 0{1) or 
smaller (where the (•) indicate a spatial average), any en- 
hancement of the photon-chameleon conversion rate due 
to electron fluctuations is expected to be sub-leading or- 
der at optical (and higher) frequencies for cluster and 
super-cluster scale magnetic fields. 



IV. CHAMELEON FIELD OPTICS 

In this section we consider how the presence of a 
chameleon field alters the properties of light propagat- 
ing through one, or many, magnetic regions. Varying the 



action Eq. ^ with respect to both 
© and 



[BF{<l>/Mo)Fn = r 



(5) 



where is the background electromagnetic 4-current: 
— 0. We consider propagation of light in an astro- 
physical background which contains a magnetic field of 
strength B. The background value of cf) is denoted 4>o{t). 
We write the perturbation in the photon field as a^^ and 
the perturbation in the chameleon field as (50. Ignoring 
terms that are 0{5<j)a^) and assuming that the proper 
frequency of the photons, w, is large compared to the 
Hubble parameter, H , we find that: 



\/S(f> X B 
M ' 
B • (V X a) 

M 



(6) 
(7) 



where 1/AI = {\nBp)^^{(j)o) and 

□00 = K=ff,0(0m,Pb,FoV4 = BV4). 

and here pb is the background density of matter. 

We assume that S(j) is small enough that we may make 
the approximation: 

1/^(00 + (50) - V4,{M = ml5^ 

where = (0o) is the chameleon mass. If the pho- 
tons are moving through a plasma with electron number 
density rie, they will behave as if they had an effective 
mass squared Wpj, where = 47rane/me is the plasma 
frequency; a is the fine structure constant and me is the 
electron mass. 



A. A Single Magnetic Domain 

We define 7_l and 7|| to be the components of the pho- 
ton field perpendicular and parallel to the magnetic field 
B, and take the photon field to be propagating in the 
z-direction. From Eq. ([6]) we have: 



■711 + 



^7_L + 



dz'^ 



B 

dzM' 
B d 



For such a system, it is well known that the probability 
of a photon, with frequency lo, converting to a chameleon 
particle (or vice versa) whilst travelling a distance L 
through a region with a homogeneous magnetic field is: 



and Af, gives Eq. where 



A 
A 

tan 26* 



Pjt^cf) — A 



sin 29 sin 



cos 26* 



m, 



off'' 



4a; ' 
2Buj 



(8) 

(9) 
(10) 
(11) 



andm^jj = m'^^-uj'^i-B^ /M'^ . Generally \B^/A'Pmlff \ < 
1 and so the last term in is dropped. Following 
Bam, we find that, up to an overall phase factor, 
7j_, 7|| and 4> = ix are transformed by passing through a 
homogeneous magnetic domain in the following way: 



7|| ^7|| 
7-L 



e'" (Vl - ^2^_L + ie-^'^ylx) 

where a = </? — A and /3 = </? + A and 

A 



tan tf — cos 29 tan 



cos 26* 



(12) 
(13) 
(14) 

(15) 



Since we must, in realistic situations, allow the light to be 
partially polarized (or even unpolarized) , it is insufficient 
to consider simply the evolution of the photon, 7j_ and 
7|| , and chameleon = 1% amplitudes. We must instead 
represent the properties of the photon field by its Stokes 
vector. We therefore make the following definitions: 



^7 

Q 

U + iV 
J + iK 
L + iM 



I. = 



(i7^n+(i7iip), 

(|7J^)-(|7||P), 
2(7±7||), 
2e''^ (711 X), 
2e''^ (7±X) • 



(16) 
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The Stokes vector for the photon field is S = 
{Ij,Q,U,V)^ , where V describes the amount of circu- 
lar polarization (CP), and Q and U describe the amount 
of linear polarization (LP). We also define the reduced 
Stokes vector: S^cd — {Q/I-y, U/I.y, V/I^)'^ . The fraction 
of light which is polarized is: 



P 



^7 ^7 

and the fractional circular polarization is: 

V 



L 



7 



We also define q — 
tion is m\ 



The fractional linear polariza- 



We normalise the photon and chameleon fiuxes so that 



A + lib 



1 (this quantity is conserved), where 



|(/)| . We also define X = 3/-y — 2. With these definitions 
we find that, on passing through a single homogeneous 
magnetic domain, the components of the Stokes vector 
transform as: 



X 



1 



-A'\X 



(17) 



-3AVi - A'^{Lsin2Lp- M cos2(p), 



(18) 



-A^/\- A^iL sin 2lp - M cos 2lp) , 
A^e-'^'iU + iV) 



(19) 



U + iV ^ vT 

Additionally, the J, K, L and M amplitudes transform 



as 



K + iJ 

L 
M 



Vl- A^e'^\K + tJ) (20) 
+Ae-'" {U + iV) , 

L cos 2(^ + M sin 2^? (21) 

(1 - 2A2)(Mcos2(^- Lsin2<^) (22) 

+aVi-a^q + x). 



From these equations it is clear that the presence of a 
light scalar field coupling to photons can result in the 
production of polarization. This is because, when a 
chameleon (or another axion-like particle) is converted 
back into a photon, that photon is polarized perpendic- 
ular to the magnetic field. If we consider the simple 
case where initially there is no chameleon flux so that 

= 1 ^ X = 1 and K = J = L = M = 0, and we 
set Q ~ Qo, U — Uo and V = Vq initially then, using 

= it is clear that upon exiting the magnetic 



domain: 



U + iV 



(1- 



i(i + 
■0)<3o 



Qo), 
1 



2^- 



P^^4.e-'"{Uo + iVo). 



If the initial total and circular polarization fractions are 
Po and Qq, their final values are 



(23) 
(24) 



where Cq = (^^(1 + Qo)V4 - A^Q^)/{1 - A^) and 
Dq = {Uq — Vq) SVC? a. ~ UqViq sin 2a. It is therefore pos- 
sible for both linearly and circularly polarized light to be 
produced. In a single magnetic domain, the production 
of the former is due to the conversion of photons into 
chameleons and then back into photons, and the latter 
is due to the birefringence of the medium which is in- 
duced by the presence of the chameleon field. If initially 
P — Pa — ^1 then after passing through a single domain: 



P = 



p2 

- P2 



We also note that if there is no initial chameleon fiux or 
polarization, = Dq = 0, no CP can be produced in a 
single magnetic domain. As we shall see below, the same 
is not true if there are multiple magnetic domains. 



B. Multiple Magnetic Domains 

In many realistic astrophysical settings, including the 
ones we will be primarily concerned with in subsequent 
sections, light passes through many magnetic domains on 
its way from a source to an observer. In each domain the 
angle, 6'„, describing the inclination of the background 
magnetic field to the direction of propagation is essen- 
tially random. Solving the full system of evolution equa- 
tions for a large number of magnetic domains involves 
diagonalising an 8 by 8 matrix as well as evaluating mul- 
tiple sums involving the random angle 0„ ^ C/[0, 27r), 
and we have been unable to find a general analytic so- 
lution, however, it is straightforward to solve the system 
numerically. This said, approximate analytical solutions 
exist in a number of interesting and important limits. A 
full presentation of the equations that must be solved in 
this set-up, and their analytic solutions in these limits is 
provided in Appendix IB] We present the results of that 
analysis below. We define N to be the number of mag- 
netic domains through which the light has passed, and in 
all cases assume that there is no initial chameleon flux. 

For fixed m'^g and magnetic domain length L, we define 
a critical frequency Wcrit such that A(a;crit) — 7i'/2, and 
hence Wcrit = When lo ^ Wcrit, P-y^^, is al- 

most independent of frequency, however when to <C i^crit j 
P-yw0 oc w^. We also define Acrit = 27r/wcrit to be the 
critical wavelength and Aosc = Adit/-^- 
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1. Weak Mixing Limit 

In a great many realistic situations we have Na <^ 1 
and NP-y^^ <C 1 and, as we shall show, the frequency 
dependence of the production of linearly and circularly 
polarized light in this limit is qualitatively similar to 
that seen in general. In this limit we must have either 
A/ cos 26* and A tan 26* < 1, or tan 26* and Atan^ 29 < 1, 
and so 

tan^ 26* , . „ , , 
a « 2A - sin 2A . 



In Appendix[B]we find that when an initially unpolarized 
light beam, with frequency w, from a single source passes 
through TV 3> 1 domains (and requiring Na <C 1 and 
NP^^^ <^ 1), the final polarization fraction, po, and 
final fractional CP, q, are essentially random variables 
and are described by the following distributions: 



NP. 



P 



7^0 



X2X4) . 



where at fixed A = m^^L/Auj, the Xi are approximately 
independent identically distributed N{Q, 1) random vari- 
ables and 



1± 



cos(2(Af- 1)A) sin27VA' 
TV sin 2 A 



When A <C ^cnt/N = Aqsc the Xi are roughly indepen- 
dent of A, but when A 3> Aqsc there is a strong A, and 
hence wavelength dependence. The above expressions 
describe the total and circular polarization fractions for 
a monochromatic light beam from a single source. If 
one has observations of many objects all at the same 
frequency the average value of p, denoted p, and r.m.s. 
average of mc, denoted q, are more useful quantities for 
comparing with observations. For the distributions above 

P = \nP,^^, (25) 



q K, \/2CT+cr_TVP, 



(26) 

where N is now the average number of magnetic regions. 
When TVA > 1, ct+ct- = 1/4 and when A^A <C 1, 
= iVA/V3. 

When some initial polarization is present (p = po 
and q = qo say, so that the initial linear polarization 
is TO/o = ^ 9o)' ^^'^ ^ different behaviour: When 
NP^^^{l-pl)/2pa > 1 (but keeping NP^^^ < 1) the 
average final polarization fractions, p and q are still given, 
to 0{Na'^, NP^^^), by Eqs. ^ and ^ respectively. 



When NP^^^{\ ~p^)/2pa ^ 1 we find instead that to 

(27) 



and 



9 = 



P = Po, 



ql{l-a^N) + 



(28) 



where this expression is only accurate to leading order in 

We shall see that for realistic astrophysical magnetic 
fields the critical wavelength Adit generally corresponds 
to UV or X-ray light. As such most polarimetry mea- 
surements of astrophysical objects will have been made 
at wavelengths ^ Aqsc = ^ait/N . For such wavelengths, 
the analytical solutions found in Appendix [B] show that 
the reduced Stokes parameters, Q/Z-y, U/I^ and V/Iry, 
exhibit a strong and oscillatory frequency dependence. 
This behaviour is very important when one wishes to 
make comparisons with observations. Polarimeters al- 
ways have some finite wavelength (A = 27r/a;) resolution, 
5\. This is usually referred to as the spectral resolution. 
Thus a measurement of the reduced Stokes parameters 
at some wavelength Aq, will actually measure an average 
of their values in the window A e (Ag — (5A/2, Aq + 5\/2). 
If one averages the reduced Stokes parameters over wave- 
length bins much larger than Aqsc much of the infor- 
mation about a chameleonic contribution will be lost. 
Specifically if we assume 5X/\ <C 1 and i5A ^ Aqsc then. 



to 0{Na , NP-y^^), p = p and rUc 



p{5y,Uj) 



Po, 
rrico 



moYi 



where: 



(29) 
(30) 



where for A^ ^ 1, Yi and I2 are independent identically 
distributed N{0, 1) random variables. Thus, in this case, 
constraints on the parameters of the scalar field theory 
could only be derived by measuring both the total polar- 
ization fraction and the CP fraction. When A ^ Aqsc, if 
the spectral resolution is too poor or the data is placed 
into too wide wavelength bins, the measured polariza- 
tion fraction, p, carries little or no information about the 
properties of (f). For optimal sensitivity to chameleonic 
effects, the spectral resolution of the polarimeter and the 
size of the wavelength bins must satisfy SX < Aqsc- We 
discuss in mV CI below how when the spectral resolution 
is sufficiently good, the strong wavelength dependence at 
wavelengths A > Aqsc can be exploited to extract strong 
constraints on chameleon-like theories from observations 
of a single object. 



2. Maximal Mixing Regime 

When A^ ^ 1 , if the chameleon to photon coupling is 
strong enough and A^A <§; 1, i.e. if A ^ Aqsc, maximal 
mixing will occur. In this limit the equations that must 
be solved simplify greatly. Further details of the calcula- 
tions are given in Appendix |Bj The strong mixing limit 
is appropriate when NP^^^ ^ 1, A^A <C 1 and A^ 1. 
When A^A <C 1 there is little production of circular po- 
larization, and so the main effect is the production of 
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linear polarization. Additionally since at high frequen- 
cies we expect qg <^ Pq for astrophysical objects we set 
the initial circular polarization fraction to zero, and the 
final CP fraction, q, remains ^ p. We find that the final 

Average Polarization Production in the Ivlaximai Ivlixing Limit 
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FIG. 1: Dependence of the mean linear polarization, p, in the 
maximal mixing limit on the intrinsic polarization (po)- The 
solid line is the exact value of p whereas the dashed line is the 
value calculated from the fitting formula Eq. p2p . The thin 
dotted line shows p = po as would be the case when chameleon 
photon mixing is weak or non-existent. We can see that for 
Po ^ 90% maximal chameleon-photon mixing increases the 
average linear polarization, whereas for po > 90% it slightly 
decreases it. 

(linear) polarization fraction, p — toi, in this limit, does 
not explicitly depend on P^^^ or any other properties of 
the chameleon field and that it is given by the following 
distribution: 

p =F{X^,cos2tP;po) (31) 
L 4(l-p2)x2 

V [(1+X2)-J,0(1-X2)COS2V^]'' 

where ip and X are independent uniform random vari- 
ables: V U[0,t:) and X ~ J7[0, 1). We note that when 
A^A ^ 1, /p does not depend on frequency. If we average 
over observations of many sources (each with A'^A <C 1 
and NP-y^fj, 3> 1) then we would measure the average po- 
larization fraction p. In the simplest case where pq = 
we have: 

1-X^ 

and so 

p = l dX^^--1.0.57. 
More generally 

p{po)^- f f dXF{X^,cos2^-po). 
I" Jo Jo 

p{po) is a monotonically increasing function of pq and 
increases from 0.57 to 1 as po goes from to 1. It is 



clear that in the strong mixing limit, the presence of a 
chameleon- like field coupling to the photon can induce a 
significant amount of linear polarization. We find that 
the following formula 

Mx(po) = ^yr^(l-2p2) (32) 
+ (|-2)(l-pg) + l. 

fits p{po) extremely well. We plot pgx against the exact 
value of p in FIG.[TJ The solid line is the exact value and 
the dashed line shows pnx- Also shown on this plot is the 
line p = Pq. For pq < 90%, the average polarization after 
maximal chameleon mixing is larger than the intrinsic 
polarization po, whereas for po > 90% it is slightly less. 
If one were to measure p < 57% for a large number of 

Probabiiity of Measuring a Given Polarization in the Maximal Mixing Limit 
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FIG. 2; Probability of measuring the linear polarization de- 
gree (p) less than some Pm for a random object if chameleon- 
photon mixing is maximal. 

astrophysical objects then at least one of the following 
must be true: NP^^^ <C 1 or A > Aosc- When A <^ Aosc, 
the chameleon induced polarization is largely indepen- 
dent of frequency, and so the spectral resolution of the 
polarimeter is not as important as it is in the weak mix- 
ing regime. The probability of measuring the total linear 
polarization of a random object to be less than some 
when A ^ Aosc and mixing is maximal (AP^^^ ^ 1), is 
shown in FIG. [H we have assumed no knowledge of the 
initial intrinsic polarization and hence marginalized over 
a uniform prior for it. 



3. General Behaviour 

When, as is often the case, one excepts little or no 
intrinsic circular polarization of the light beam i.e. mco = 
0, we are able to combine the results presented in the 
previous two subsections to provide a fitting formula for 
the general form of the mean value of p after the light 
beam has passed through A 1 magnetic domains. We 
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find 



p{N) 



where 



\ 



In the maximal mixing limit p = p^^{po). In the weak 
mixing limit, when NP^^^ <^ 1, we have: 



piN) 




bHp^HPo)~pl)N^P^^^ 



So if Pq is small, we have p{N) w as required. 

If instead po is larger, p{N) = _po + 0{N'^P^^^). This 
provides a very good fit to the simulated data in all cases. 



C. Optical Signatures of Chameleon Fields 

We presented above the results of a mathematical anal- 
ysis of how the presence of a light scalar field coupling 
to matter would alter the polarization properties of light 
passing through a magnetic field, the details of which 
can be found in Appendix [B) By combining the results 
of this analysis with numerical simulations, we now de- 
tail the main physical signatures that a chameleon field 
would imprint on the polarization properties of light from 
astrophysical sources. Above we found that there were 
two main effects: 

• the production of polarization, 

• the production of circular polarization. 

Each of these effects depends on frequency in a charac- 
teristic manner that is well illustrated by considering the 
weak mixing limit of ijlVB II above with no initial polar- 
ization (po = 0). In this hmit > 1 but NP^^^ < 1 
for aU A and Na^ < 1. This requires BL/2M < 1; 
f'max = limA^o-P7^0 ~ {BL/2MY <C 1 is the maxi- 
mum value of P-y^0. 

We consider the weak mixing limit, assuming that 
the polarimeter has wavelength resolution < Aosc = 



Acrit/A^ = 47rV|TO^jj|Lpath where 



'eff 



1^ - Wpi 



and 



-f-path is the total path length of the light through the 
magnetic field. In this limit, when there is no initial po- 
larization, both the induced degree of polarization: p, 
and circular polarization q = \mc\ are proportional to 
NPj^^. FIG. [3] shows possible simulated forms for the 
rescaled total polarization degree, p/NPamx, linear po- 
larization TOi/A^-Pmax and CP q/NP^s,x for two different 
hypothetical objects, where for example A^ « 100 in both 
cases. We can clearly see from this that production of lin- 
ear polarization is greatest for A < Acrit and CP polariza- 
tion production is peaked in the region Aosc ^ ^ ^ Adit- 




10 10 10 10 

Wavelength:X/X 




FIG. 3: Dependence of the total polarization degree, p, the 
linear polarization degree, mi and the circular polarization 
degree q on wavelength for two hypothetical objects with A'' = 
100 and A/'Pmax < 1. Here Acrit = 47r^/|m^fl|L where L is 
the coherence length of the magnetic field regions. The total 
path length of the light through the magnetic field is given 
by I/path = NL. We define Aosc ~ Xait/N. We have assumed 
that initially p = and that there is no initial chameleon flux. 



As expected, we can also see that both the induced lin- 
ear and circular polarization degrees are highly frequency 
dependent for A > Aosc- 

Assuming 6X < Aqsc, averaging p, and q over many 
sources each at roughly the same distance, and hence 
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FIG. 4: Dependence of the total average polarization de- 
gree, p, and r.m.s. circular polarization degree, q, when av- 
eraged over many sources, each with A'^ = 100 and A'^Pmax ~ 
N{BL/2Mf < 1. Here Acrit = 47r^/|m^fj|L where L is the 
coherence length of the magnetic field regions. The total path 
length of the light through the magnetic field is given by 
ipath = NL. We define Aosc = Xait/N. Initially we have 
assumed that p — and that there is no initial chameleon 
flux. 

with roughly the same N, gives p and q. The forms of 
p, and q are shown in FIG. HI We see that both quan- 
tities grow strongly when A ~ Acrit and that q is peaked 
between Aosc and Acrit- The form of the averaged CP 
degree, q, is very distinctive. The height of the peak, 
as well as the maximum value of p are determined by 
N{BL/2M)'^, whereas the position of the peak and its 
width are fixed by Aosc and Acrit- If such a peak should 
be resolved, one would in principle be able to determine 
both N{BL/2MY , TTioff and the coherence length, L, of 
the magnetic field regions. Measurements of circular po- 
larization for Aosc J; -^crit could therefore provide 
a powerful tool with which to constrain chameleon theo- 
ries; we discuss this further in Will Qualitatively similar 
behaviour is seen for the chameleonic induced polariza- 
tions in the more general case where iVPrnax can take any 
value and the intrinsic (i.e. non-chameleonic) polariza- 
tion is not restricted to vanish. When A^Pniax ^ 1 is 
allowed, the chameleonic production of linear polariza- 
tion is peaked for A < Amax where: 

Amax BL \ 

= max 1, — — — . (33) 

Acrit V ttVNMJ 

CP production is in general peaked between Aosc and 

Amax- 

In practice, however, as we shall discuss further in WII 
below, it is rare for current polarimeters to have 5X <C 
Aosc; although measurements do exist with 5\ ~ ©(Aqsc)- 
As well as requiring 5\ < Aosc , to measure p, rhi and q one 
must also have measurements of many sources where the 
light from each source is expected to have passed through 
roughly the same number of magnetic regions, each with 
roughly the same properties, as the light from any other 
source. This requirement introduces a fair amount of un- 



certainty and will ultimately limit ones ability to accu- 
rately constrain the averaged quantities. Another prob- 
lem is that even when the intrinsic polarization is small 
(po ^ 1), if NP^ <C 1, the form of p, and q are highly 
dependent on pq. Unless one can measure or accurately 
predict the intrinsic polarization, this again limits ones 
ability to accurately constrain chameleon theories. 

Many astrophysical polarization measurements are 
made at wavelengths A ^ Aosc where the chameleon in- 
duced contribution to the Stokes parameters exhibits a 
highly oscillatory wavelength dependence. When po 'ti 1 
and provided SX ~ O(Aosc) or smaller, we can exploit 
this property to extract strong constraints about the 
properties of any chameleon-photon interaction from ob- 
servations of a single object without any detailed prior 
knowledge of pq. We do this by defining a smooth- 
ing scale (5Asmooth which is picked to be ^ Aosc but 
smaller than the wavelength scale over which the in- 
trinsic polarizations, po, mio and qo, are expected to 
vary strongly. By removing the smoothing signal from 
the measured signal we should recover a superposition of 
any induced chameleonic signal and the noise. Assuming 
that the noise is either random or that it does not have 
a wavelength structure that mimics that of the induced 
chameleon signal we can extract constraints on M. Fur- 
ther details of how M can be constrained in this manner 
are given in Appendix [Cl Using this method, it is possi- 
ble to extract strong constraints on M using data from 
only a single source. 



V. LARGE SCALE ASTROPHYSICAL 
MAGNETIC FIELDS 

The largest scale magnetic fields that are known to ex- 
ist are those associated with galaxies and galaxy clusters. 
In both cases the mean field strength has been measured 
to be roughly a few micro Gauss. It is also thought likely 
that a weak, B < 10~^G, magnetic field permeates the 
inter-galactic medium (IGM). We discuss the observed 
properties of first two fields as well as the hypothesised 
properties of the latter below. The electron density, rie, 
determines the plasma frequency, ujpi, which plays a crit- 
ical role in determining the effective chameleon mass, 
TO^g, and hence also critical frequency, Wcrit, above which 
polarization production is peaked. Therefore we also 
quote the observed or estimated values of Ug for each 
of the three regions. 



A. Galactic magnetic fields 

Galactic magnetic fields, particularly those of our own 
galaxy, could produce detectable polarization effects if 
chameleon-like fields interact strongly enough with pho- 
tons. Galactic magnetic fields have been observed to be 
a superposition of a regular magnetic field, Breg, and 
a random magnetic field. Brand (see [sJl and references 
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Brcg — 2 fiG with field reversals occurring over kilopar- 
sec scales, and -Brand = 3^G with Lrand — 20 pc. The 
average electron density was taken by Sun et al. to 
be ric = 0.03 cm~'^. Generally light observed from ob- 
jects within our own galaxy will have passed through 
N ^ 0(1) regions of the regular magnetic field, but 
N ^ 1 different coherent regions of the random mag- 
netic field. Taking L = ij-and = 20 pc and B = 3fiG for 
the random magnetic field and Uc — 0.03 cm^'' we have: 



and 



therein). The regular component of the magnetic field L 
has a coherence length, I/^eg ~ few kpc i.e. about the 
scale of the galaxy [5l|. The component of the regular 
part of the magnetic field along the line of sight to distant 
objects such as pulsars and extragalactic radio sources 
has been measured using Faraday rotation. These mea- 
surements are performed using electromagnetic waves 
whose frequency is well below Wcrit- The interpretation 
of such measurements would therefore be largely unal- 
tered by the presence of a chameleon field or similar 
light scalar field. The average regular magnetic field in 
own galaxy is locally (within about 2 kpc of the Sun): 
Brog = 1.8 ± 0.4 ^G gail], rising to about 4.4±0.9/iG 
in the more central Norma arm [54j . The magnetic field 
is aligned with the disk of the galaxy, and is coherent out 
to a galactic radius of about 5 kpc, field reversals then 
occur at i? = 5 kpc, 6 kpc and 7.5 kpc [Sll. [ssj. 

The random magnetic field. Brand, is often slightly 
larger than the regular magnetic field. The largest scale 
of the turbulent field was determined from pulsar rotation 
measures (RMs) as irand = 55 pc by Rand and Kulkarni 
[5^, with a turbulent field strength about about 5/xG. 
A similar study by Ohno and Shibata [STj found Lrand = 
10 100 pc with a random field strength of 4 — 6 fxG. Lrand 
has also been estimated by the depolarization of light by 
turbulent fields at centimetre radio wavelengths, an d by 
Faraday dispersion at decimetre radio wavelengths [55|. 
Both methods give results consistent with Lj-and ~ 20 pc. 

Recently Sun et al. combined radio tele- 

scope and WMAP measurements of diffuse polarized ra- 
dio emission from the Milky Way with Faraday rota- 
tion measurements to obtain an overall model of the 
Milky Way's magnetic field. They found that on average 



2 kpc we have: 
\B\L 



2M 



0.612 



10^° GeV 
M 



,(rcg) 



VffI 



27r 



= 2.04 keV 



(36) 



(37) 



so when | 
note that 



_ , ,2 \(rcg) 
— ^plJ -^crit 

< cjHi for nrirh 



= 27r/wcrit = 6.08 A. We 
< 1.3 X 10-" eV. 



VffI ""pi 

In the weak mixing regime, i.e. when the chameleon in- 
duced polarization is small, we find that the total induced 
polarization is a sum of that which would be separately 
induced by the random and regular magnetic fields. 

For the dark energy inspired chameleon model dis- 
cussed in i]II Al we have < v^Wpi in the galaxy, and 
hence \ml^\ < w^i, for aU Mq > 3.9 x 10^ GeV when 
n < 3.3. 



B. Intracluster magnetic fields 

In galaxy clusters electron densities of ric w 10^^ cm"^ 
are typical, as are magnetic field strengths of a few fiG, 
rising to tens of /iG at the centre of cooling core clusters. 
These magnetic fields are coherent over length scales of 
about L fa 10 ~ 100 kpc [11]. Galaxy clusters typically 
extend over a length scale of idust ^ 1 Mpc. A light 
beam traversing a galaxy cluster would therefore pass 
through roughly N = Lc\ust/L w 100 — 1000 magnetic 
regions. 

In a study of data from 53 radio sources located in and 
behind Abell clusters and a control sample of 99 sources 
Kim et al. 59] found the mean core electron density of 
a cluster to be ric — 3.5 ± 2.7 x IQ-'^cm"'^, where the 
radius of the core is rcoro = 0.65 ± 0.41/i~^Mpc, and 
found cluster magnetic field strengths of 0{1)^G with 
coherence length ~ 10 kpc. A study of 18 radio sources 
close in angular position to the Coma Cluster by Kim et 
al. [63] found the following result for the strength of 
magnetic fields in the intracluster medium (ICM): 



^) = 0.92 X 10- f 1^^) , (34) 
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1/2 
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10 kpc 
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(38) 



and Wpi = 6.4 x 10 eV so 



,(rand) _ l"^gff]L 



27r 



= 20.4 eV 



'off I 



Pl 



When <^ 6.4 x 10 eV and hence 



, (rand) „ / 

A„.,^ = 27r/a;crit 



608 



ffl - '^pl' 

— i^cm — For an object in our galaxy 

at a distance d we take N k, (i/20pc. Therefore if, as is 
typical, d ^ \ kpc we have N k, bO. 

Taking B = 2 fiG for the regular magnetic field and 



where /175 is defined in terms of the Hubble parameter 
today: Hq = 75/175 km s~^ Mpc^^. A subsequent study. 



again of the Coma cluster, by Feretti et al. [61| found 



(35) tangled magnetic fields with length scales of about 1 kpc. 



1/2 

so i?icM ~ 7.9/175 . A study of 16 low redshift (z < 
0.1) "normal" galaxy clusters by Clarke, Kronberg and 
Bohringer [g^] found that the ICM of these clusters was 
permeated by a slightly larger magnetic field: 
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Based on the studies of Kim et al. [59|, |60| , we take 
the foUow representative values for the parameters which 
describe magnetic fields in the ICM: 



and 



1 kpc, 



B ^ 7.9h 



1/2 
75 ' 



3.5 X 10^3 cm~^ 



We define Lpath to be the path length a given light beam 
traverses through a cluster, and take as a representative 
value Lciust = 1 Mpc. The number of magnetic regions, 
A'', is given by = ^path/^, and hence we take 

N = 1000. 

With these values we have ujp\ = 2.2 x 10^^^ eV and 



\B\L\ 
2A/ JjcM 



and 



(ICM) 



2tt 



= 120cV J— p . (40) 



-'pi 



When \mlff \ = cj^j, A^f^^' = 27r/cjcrit = 104 A. We note 
that < when < 4.4 x 10^^^ eV. 

In galaxy clusters, we have < V^i^pi for the 
chameleon model introduced in mi Al for all > 3.9 x 
lO*' GeV when n < 3.5. 



C. Intergalactic magnetic fields 

Although a number of different mechanisms have been 
suggested that would produce large scale magnetic fields 
in the intergalactic medium (IGM), at the present time 
very little is known about whether such fields actu- 
ally exist, let alone their typical strengths. A coherent 
magnetic field on the current horizon scale would pro- 
duce an anisotropic expansion. CMB and Faraday ro- 
tation constraints on such a scenario limit B < 10~^ G 
[63l . H^l . Faraday rotation also constrains smaller scale 
magnetic fields. For a 50 Mpc coherence length one has 
B < 6x 10-9 G, and B < 10"^ G for Mpc scale co- 
herence lengths The CMB has also been shown to 
constrain fields with a coherence length between 400 pc 
and 0.6 Mpc to be < 3 x 10"^ G Motivated by 

the need to explain the origin of galactic magnetic fields 
it is thought that IGM magnetic fields with coherence 
lengths of a few Mpc are likely (see and references 
therein). Most of the proposed theoretical mechanisms 
for generating such fields would, however, only produce 
them with strengths well below the current observational 
upper bounds [6Q|. These seed fields are then amphfied 
by some dynamo mechanism during galaxy formation to 
the ~ ^G galactic magnetic fields observed. 

Typical electron densities in the IGM are Ue ~ 2.5 x 
10^''cm~'^ giving Wpi = 1.8 x 10~^'*eV and so 



B 



BL /10i''GcV\ 
2M - '■''\-^)\lO-^G 



L 



IMpc 



(41) 



Wcrit = 3.4 eV 




hence for 



'off I 



and L — 1 Mpc, we have A, 



■crit 



3647A. For the dark energy chameleon potentials dis- 
cussed in mi Al the mass of the chameleon due to the den- 
sity of the IGM is < \/2[jpi i.e. I^Ti^jjl < tjpj for n < 4.5 if 
Mq > 3.9 X 10^ GeV and the chameleon couples only to 
baryons. If the chameleon couples to dark matter with 
equal strength then the same is true but only for n < 3.5. 



VI. CURRENT POLARIZATION 
CONSTRAINTS ON CHAMELEON-LIKE 
MODELS 

In this Section we review a number of astronomical 
polarization observations and deduce how they constrain 
the properties of any chameleon-like field. 

We noted in ijlVB II that at wavelengths A < Aosc = 
Acrit/A^j any chameleon induced polarization signal is a 
highly oscillatory function of wavelength, with oscilla- 
tion length « Aosc- This is particularly important as 
many astrophysical polarization measurements are made 
at optical frequencies for which A < Aosc, and either 
the Stokes parameters are put into wavelength bins with 
width 3> Aosc , or the spectral resolution of the polarime- 
ter is so poor that it effectively averages over a range of 
wavelengths which is 3> Aosc- In either situation, any sig- 
nal of chameleon mixing will be washed out, and no con- 
straints on the chameleon model are possible. If there is 
no initial polarization the polarization fraction depends 
on A via Pj^^ and cr^., these are all highly oscillatory 
functions of A except when A < Aosc- We note that in 
some cases the spectral resolution of the polarimeter is 
good enough to resolve any chameleon induced polar- 
ization, but the published data only quotes the Stokes 
parameters in bins much wider than 6X. In these cases, 
the published data cannot bound chameleon-like theories 
but constraints should follow from a reanalysis of the raw 
data. By way of an example, we have performed such a 
reanalysis for observations of three stars in our galaxy, 
however in general such a reanalysis is beyond the scope 
of this article, and is intended to form the basis of a fu- 
ture work. 



A. Starlight Polarization 

Polarization is not usually produced by the thermal 
emission of stars. In [gJ a statistical analysis of the 
larg est available compilation of galactic starlight data 
[68| was performed. The data is statistically significant 
for sources out to distances of 6 kpc, and the average po- 
larization of light from stars at such distances is 2%. 
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This data, provided in the polarization catalogue [68l |. 
is in very wide wavelength bins that generally cover the 
whole range of optical frequencies, i.e. the bin width 
is SX ~ 2000 — 8OOOA. For comparison, the oscillation 
length, Aosc — Xcrit/N, in the galaxy for such stars is 
Aosc ~ 2 — 12 A for both the random and regular com- 
ponents of the magnetic field. Thus SX 3> Aosc and the 
data provided in [68] as well as the subsequent analysis of 
[67j does not provide useful constraints on chameleon-like 
theories. 

Existing starlight polarization measurements can, how- 
ever, constrain chameleon-like theories. UV polar- 
ization of starlight was measured for 121 objects by 
the Wisconsin Ultraviolet Photo-Polarimeter Experiment 
(WUPPE), which flew on the ASTRO-1 and ASTRO-2 
NASA space shuttle missions, and had a nominal spec- 
tral resolution of 6A[6^. The data from these observa- 
tions is available from the Multimission Archive at STScI 
(MAST) [70]. A full rcanalysis of the data for aU 121 ob- 
jects is beyond the scope of this work, however we have 
derived preliminary confidence limits on BL/2M in the 
galaxy using data from three objects: HD2905, HD37903 
and HD34078. These objects were picked as they all lie 
at distances between 500 pc and 1000 pc, which is not 
so close that a chameleonic signal would be too small to 
detect, and not so far away that numerical calculations 
involved in extracting the confidence limits on BL/2M 
are too time consuming. Other than that, the choice of 
objects is entirely arbitrary. A detailed account of the 
method and resulting confidence limits derived from the 
polarization measurements of these objects is given in 
Appendix [Cl Intriguingly, we found that the data from 
all three objects preferred a non-zero value of BL/2M at 
a prima facie statistically significant level. This analy- 
sis shows that there is some structure in the polarization 
data which is consistent with the signal that we predict 
would be induced by a chameleon field. It would be pre- 
mature, however, to claim this as an actual detection 
before a similar analysis has been conducted for more ob- 
jects, and before a thorough analysis of all systematics 
which could be sources this signal has been undertaken. 
The most conservative, in the sense that they are the 
widest and are expected to be the most robust, confi- 
dence intervals were found using the bootstrap-t method 
(see Appendix [C] for further details). At 95% confidence 
we found, taking i w 20 pc and Acrit = 



mean 6.27 x 10 ^ and variance cr^; a = 0.58 x 10 ^. 
Hence we find the following approximate confidences 



\B\L 
2M 
\B\L 
2M 
\B\L 
2M 



rand 



(4.68ti;7o) X 10^^ (HD2905), (42) 
(7.59ti:^^) X IQ-^ (HD37903), (43) 
(8.58j:?:^^) X 10"2(HD34078).(44) 



If we assume that the same value of \B\L/2M is ap- 
propriate for each object, by combining the polariza- 
tion data for all three stars we find that the estimate 
of BL/2M is approximately normally distributed with 



\B\L 
2M 

\B\L 

2M 



= (6.27 ± 1.14) X 10"^ (95%), (45) 



rand 



= (6.27 ± 1.91) X 10^2 (99.9%). (45) 



From this preliminary analysis, it therefore appears as if 
the polarization data of the three objects considered is 
consistent with a value of BL/2M which deviates from 
by more than lOcr. 

Although this analysis is only preliminary, it does 
appear as if there is a reasonably significant, and ro- 
bust, statistical preference towards the existence of a 
chameleon-like field in the starlight polarization data of 
the three objects we have considered here. This is a 
highly surprising result, and as such it would be pre- 
mature to claim it as a detection. Whilst it is well 
beyond the scope of this particular article, a thorough 
analysis of possible backgrounds and sources of system- 
atic error which could mimic the signal from a chameleon 
field would have to be undertaken before any such claim 
could be made with true confidence. In particular, since 
all the data analysed comes from a single experiment 
(WUPPE) it possible that the 'detection' of a non-zero 
value for BL/2M is actually due to effects intrinsic to the 
instrument. In order to quantify the magnitude of such 
instrumental effects it would be necessary to study sim- 
ilar data from other polarimeters. We have only consid- 
ered three of the well over one hundred objects measured 
by WUPPE. When more objects have been analyzed it 
should be possible to better estimate the effect of sys- 
tematic error in the determination of BL/2M by consid- 
ering the spread in the values of BL /2M determined for 
each object. What we can say with confidence is that 
there is some structure in the polarization of three ob- 
jects considered which is not consistent with either ran- 
dom error or that predicted to be induced by interstel- 
lar dust. Furthermore this structure exhibits non-trivial 
oscillatory frequency correlations which at least in part 
mimic that predicted by the chameleon model. At the 
present time we cannot rule out possible systematic ef- 
fects having a relative magnitude of 0(1). The presence 
of such effects would raise both the extracted upper and 
lower bounds bounds on BL/2M. Whilst this means 
that any non-zero lower bound on BL/2M can only be 
seen as tentative at best, the upper bounds on BL/2M 
should be robust. We therefore believe it to be better to 
see the data as providing the following 95% and 99.9% 
confidence upper bounds on BL/2M: 
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< 7.2 X 10^2 (95%)^ (47) 

1 

< 8.1 X 10^2 (99.9%). (48) 



We also consider observations of the UV polarization 
of two stars made with the Faint Object Spectrograph 
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(FOS) of the Hubble Space Telescope (HST) and re- 
ported in Ref. [7l|. In this case, we have only un- 
dertaken a preliminary analysis of the data, postpon- 
ing a full reanalysis to a later work. Observations were 
made for 1279A < A < 3300A, and the HST FOS has 
a nominal spectral resolution of 2 — 4A in this range. 
The data published in Ref. fzi'l was binned to give ten 
data points in each frequency region they considered, al- 
though the precise width of the bins is not stated. The 
shortest wavelength region was 1279 - 1603 A. Assum- 
ing that each of the ten bins in this region had equal 
width, the bin width is 32.4A. The two stars, HD7252 
and HD161056, are respectively 824 pc and 295 pc from 
the earth. This gives an oscillation wavelength, Aosc = 
Acrit/-/V, no smaller than 15A and 41 A for the random 
and regular components of the galactic magnetic field 
respectively, when < 9 x 10^^^ eV. A significant 
amount of a chameleonic signal should therefore survive 
the rebinning process in the 1279 — 1603A wavelength 
grating; this may not be the case for the lower frequency 
gratings. In the 1279 — 1603A grating, the polarization 
angle of HD7252 was found to be independent of fre- 
quency with a standard deviation of about 5 degrees. It 
is noted in Ref. [TlJ however that the systematic uncer- 
tainty in the polarization angle could be 10 degrees or so. 
This corresponds to the component of the reduce Stokes 
vector, P± say, that is perpendicular to the mean polar- 
ization detection in the region 1279 - 1603A satisfying 
|Pl| < 0.2%. Assuming < 9 x 10~^^ eV in the galaxy 
and that Lrand ~ 0(20 pc) (the precise value of Lrand 
does not greatly alter the resulting constraint) so that 
Acrit = 608A, and using the method outlined Appendix 
IC 2[ we find the following 95% and 99% confidence limits 



BL \ 

— < 8.9 X 10^2 (95%), 
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\2M) 



< 12.7 X 10^2(99 9%)^ 
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These constraints are consistent with, but weaker than, 
those found from the WUPPE data. 



B. The Crab Nebula 

The polarization of X-ray light from the Crab nebular 
was reported in [72]. The measured linear polarization 
fraction was p = 18 ± 4% at a frequency of = 5.2 keV 
and p = 16±2% at a frequency of 2.6 keV. This confirmed 
the hypothesis of synchrotron X-ray emission. The Crab 
nebula is at a distance of 2kpc from the solar system 
so photons from the Crab nebular pass through 0(1) 
regular magnetic domains and 0(100) random magnetic 



domains to reach the earth. When 



= 6.4 X 



10-12 gv, we have Aosc « 6A; iVwcrit ~ 2 keV. Thus both 
measurements are in the u> > iVwcrit region, where p is 
almost independent of frequency. The spectral resolution 
of these measurements is SX < 0.7A<C Xcrit/N. 



The linear polarization fraction, p is given by a prob- 
ability distribution even if we are in the maximal mixing 
regime and so the amount of information one can extract 
from a single measurement is limited. We found that the 
average polarization fraction for a set of objects in the 
maximal mixing limit is > 0.57. However, if po <^ 0.16 
initially, one would still expect to measure p < 0.16 — 0.18 
for a given object about 17% of the time. Even the pos- 
sibility of maximal mixing at X-ray frequencies cannot 
therefore be ruled out by the Crab Nebula data. 



C. Type la super novae 

In [tI, [tII supernova polarimetry data published be- 
fore 1996 was studied. The degree of polarization of light 
from type la supernovae was less than 0.2 — 0.3%. In 
[t^ high-quality spectro-polarimetry data was reported 
for SNIa 2001cl. It was found that the maximum lin- 
ear polarization of the light from the supernovae was 
p w 0.2 — 0.3%. The supernova was observed at frequen- 
cies w « 1.4 eV— 3.8 eV and the spectral resolution of the 
polarimeter was 5X w 12. 7A. The Stokes parameters were 
later re-binned into 5X = 15Abins. The supernova lies at 
a redshift of z w 5 x 10"'^ corresponding to a distance of 
roughly 20 Mpc. 

If, as light travels from the supernova to the earth, 
mixing with chameleons occurs mostly in the intergalac- 
tic medium (as opposed to in galaxies or clusters) then 
the PVLAS bound rules out maximal mixing. The 
critical frequency for weak mixing in the intergalactic 



medium is lu, 



(IGM) 



3.4(Ligm/ Mpc) eV, where L is 



the coherence length of the IGM magnetic field. Since 



NLiGM = 20 Mpc we have Nuj^ 



IGM 



68 eV; Ao 



182A. Hence SNIa 2001el was observed at frequencies 
u! <^ NluI^^. Any chameleon induced polarization frac- 
tion would therefore be a highly oscillatory function of 
the frequency. For this chameleon signal to survive 
the binning process and be detected one must ensure 
that the polarimeter's spectral resolution and width of 
the wavelength bins satisfy 6X < 182A. In this case 
SX — 15Aand so the data can indeed be used to con- 
strain chameleon- like theories. At the wavelengths ob- 
served, the chameleonic signal would look like random 
noise that grows with frequency. A full analysis of the 
raw data reported in [t^ is beyond the scope of this 
work. However, a preliminary analysis of the scatter 
in the component of the Stokes vector perpendicular to 
the mean direction of polarization, P±, in the frequency 
range 4181 — 863lA provides strong constraints. At five 
different epochs, it was found that P± was consistent 
with zero to within about ±0.3%. We take \Pj_\ < 0.3% 
and extract approximate upper confidence limits on the 
chameleon to photon coupling using the method detailed 
in Appendix IC 21 When < 2.5 x lO-^^eV in the 
IGM, and Ligm O(lMpc), we find the following 95% 
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and 99.9% confidence limits 



< 5.2 X 10-^ (95%), 



IGM 



\B\L 
2M 

< 7.2 X 10-2(99.9%). 

/ IGM 

If the intergalactic magnetic field is sufficiently small 
(i.e. B < IQ-^^ G) mixing between light from the super- 
nova and chameleons will occur mostly in galaxies and 
galaxy clusters. SNIa 2001el is located in the nearly 
edge-on spiral galaxy NGC 1448, however the line of sight 
does not intersect with either the core or the disk of the 



host galaxy [75| . Additionally, at only 20 Mpc away, light 
from SNIa 2001el does not pass through any significant 
intra-cluster magnetic fields. Our solar system currently 
lies close to the midpoint of the galactic plane, and mod- 
els of the galactic magnetic field and electron density sug- 
gest that it has a scale height above the midpoint galactic 
plane of about a kiloparsec. At the very least then, light 
from SNIa2001 will have passed through roughly 1 kpc 
of the random galactic magnetic field. For the random 
galactic magnetic field we have Xl^^t'^/N « 12 A, where 



N 



50 when |m 



effi 



pi 



6.4 X 10-12 eV. There- 



fore SX ~ AJJ^j"'^/iV « I2A, and a chameleon signal could 



be detected. Since SX 



K7^/^^ the chameleon in- 



duced polarization would look like random noise. Again 
a preliminary analysis of the data of [75] , gives the follow 
95% and 99% confidence limit, where we have assumed 
< v^Wpi « 9 X 10-12 eV in the galaxy and that the 
coherence length of the random component of the galaxy 
magnetic field is 0(20 pc) so that Adit ~ 608A: 
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< 0.14(95%), 

< 0.18(99.9%). 



D. High Redshift Quasars 

The optical and UV polarization of some high red- 
shift quasars have been measured [ll, [t^] often using 
the HST FOS. Below frequencies w ~ 1 eV the quasars 
have a polarization of about 1% but there is an interest- 
ing rise in the polarization above frequencies w « 2.5 eV. 
At electron- Volt frequencies mixing between photons and 
chameleons is expected to be highly frequency dependent. 
The HST FOS has a nominal spectral resolution of 2— 4A, 
which is in principle good enough to resolve the expected 
chameleon signal if rricf, <C 6.4 x lO-^^ eV in galaxies or 
galaxy clusters. The data in Jt^-Itt!] is then rebinned with 
bin widths of SX — 32 — 270A. Extracting the most strin- 
gent constraints on chameleon theories would require a 
full reanalysis of original data. This is beyond the scope 
of this article. However, by analysing the data of Impey 
et al. [zil for object PG 1222-1-228 at z « 2 as rebinned 



and presented in Ref. [77[ we can extract useful con- 
straints. Specifically, we focus on the spread of the Stokes 
parameter that is perpendicular to the mean polarization 
angle. The light from this QSO will have travelled at least 
« 1 kpc. We assume that the coherence length of the ran- 
dom component of our galaxy's magnetic field, Liand, is 
0(20 kpc). Making the conservative assumption that the 
total path length through our galaxy's magnetic field is 
Ikpc, when < V2ijJpi ~ 9 x lO-^^eV in the galaxy, 
we find the following 95% and 99% confidence limits: 



\B\L 
2M 
\B\L 
2M 



< 0.6(95%), 

< 1.1(99.9%). 



We expect that a full reanalysis of the original data would 
raise this limit greatly as currently the bounds are consid- 
erably weakened by the relatively large size of the wave- 
length bins (compared to Aqsc)- 

If there is a sufficiently strong intergalactic magnetic 
field then this would also produce chameleon-photon mix- 
ing. We make the conservative assumption that IGM 
magnetic fields only go out as far as z = 1, so that the 
propagation distance through the IGM magnetic field 
is about 2.5 Gpc. Assuming the IGM magnetic field 
is coherence over roughly megaparsec scales, we find 
the following confidence limits when < v^^pi ~ 
2.5 X 10-1'' eV in the IGM: 



\B\L 
2M 
\B\L 
2M 



IGM 



IGM 



< 1.4 X 10-2 (95%), 



< 2.1 X 10-2 (99.9%). 



These constraints are particularly strong because the 
quasar is so far away, and as such the light from it travels 
through many different coherent regions, ~ 0(2500), of 
any IGM magnetic field. This counter balances the loss 
of information due to the relatively large width of the 
wavelength bins. A full reanalysis of the raw data would 
likely raise these bounds on M. 



E. Gamma Ray Bursts 

Measurements of linearly polarized gamma rays have 
been made for four GRBs and these are summarised in 
Table HI The last observation has been challenged ^] . 

TABLE I; GRB Polarization Measurements 



GRB930131 


[79] 


0.35 <p< 


1 


3 keV < uj 


< 100 keV 


GRB960924 


[79] 


0.5 <p < 


1 


3 keV < UJ 


< 100 keV 


GRB041219a 


[80] 


0.56 <p< 


1 


100 keV <Lu < 350 keV 


GRB021206 


m 


0.6 <p < 


1 


0.15 MeV < 


UJ <2 MeV 



GRBs are the only objects we consider that are believed 



18 



to be highly polarized initially. Theory predicts the emis- 
sion of highly linearly polarized light with 0.6 < toi < 0.8 
due to synchrotron emission. This hypothesis was con- 
firmed by observations of polarization in the GRB after- 
glow [8^. 

Mixing at gamma ray frequencies is maximal in all 
galaxies and clusters regions if M ^ 10^ GeV: maximal 
in galaxies if M < few x 10^ GeV and in the ICM if 
Af < 4 X 10" GeV. li B ^ IQ-^ G in the IGM, then 
maximal mixing would occur if Af < 5 x 10^° GeV; how- 
ever if TO^ <C 2.2 X 10~^^ eV in the IGM then this scenario 
is strongly ruled out by the bounds obtained above. If 
the mixing is maximal the mean observed linear polar- 
ization at high frequencies should he p > 0.57; consis- 
tent with current GRB observations. It is not possible to 
make a more precise prediction than this without know- 
ing more accurately the initial polarization of the GRB. 
A better understanding of the central engine of the GRB 
and better polarimetry for GRBs would allow strong con- 
straints to be placed on the chameleon model. If future 
observations constrain p < 0.57 maximal mixing in the 
chameleon model would be ruled out, and strong con- 
straints on M would follow. If p > 0.8 is observed such 
a high degree of polarization cannot be explained by the 
synchrotron mechanism and a chameleonic explanation 
would be favoured. 

If M is very large the mixing between chameleons and 
light from GRBs would be weak. Then it becomes diffi- 
cult to put bounds on the chameleon model both because 
of the limitations of polarimeters and, if there is no inter- 
galactic magnetic field, the difficulty of estimating how 
many magnetic domains have been traversed. 



F. CMB Polarization 

The upcoming Planck satellite will measure the polar- 
ization of the CMB to a high degree of accuracy. How- 
ever it is extremely hard to estimate how many magnetic 
domains radiation from the CMB would have passed 
through, particularly as if there is an intergalactic mag- 
netic field it is not known whether this field is primor- 
dial. Neglecting the intergalactic magnetic field it might 
be possible to use galaxy and cluster surveys to estimate 
how many magnetic domains the radiation had passed 
through. However because the frequency of CMB radia- 
tion is so low mixing with the chameleons will be weak 
and highly oscillatory and the amplitude of these oscil- 
lations is damped as w^. A weak and highly oscillatory 
chameleon signal would be very hard to detect. 



G. Summarised Constraints 

The tightest constraints on the chameleon to matter 
coupling come from the WUPPE starlight polarization 
data, in the context of photon to chameleon conversion 
in the galaxy, and from HST FOS measurements of the 



polarization of high redshift quasars in the context of 
conversion in the intergalactic medium. Our preliminary 
analysis of starlight polarization data appears to provide 
a non-zero lower bound on 1 /Af, however for the purposes 
of this discussion we only consider the upper bounds on 
1/Af here. For the IGM we took the coherence length, 
L, to be 1 Mpc and for the galaxy we assumed, L = 
20 pc, however the precise values of these quantities do 
not greatly effect the upper bounds on BL/2M . Taking 
these typical values for L and B « 3 /xG for the strength 
of the random component of the galactic magnetic field, 
we find at 95% confidence: 



Af > 1.3 X 10^ GeV, 



Af > 1.1 X lO^GcV 



B 



IGM 



10-9 G 



At 99.9% confidence we find similarly 
Af > 1.1 X 10^ GeV, 



Af > 7.3 X 10^° GeV 



10-9 G 



(49) 
(50) 

(51) 
(52) 



In both cases the upper constraint applies if < 
1.3 X IQ-^^ eV in the galaxy and the lower one if < 
2.5 X IQ-^^eV in the IGM. Since i?iGM is currently un- 
measured, the strongest constraint on Af is comes from 
the starlight polarization measurements, the interpreta- 
tion of which relies only on knowledge of the galactic 
magnetic field. If i^icM ^ IQ-^^ G, however, then the 
constraints coming from high redshift quasars currently 
provide the tightest lower bounds on Af. In terms of 
chameleon theories, these constraints represent an im- 
provement of almost 2.5 order of magnitude on the pre- 
vious best lower bounds on M coming from laboratory 
tests, specifically M > 3.9 x 10^ GeV at 99.9% confidence 
from the GammeV experiment [2]| . 

GammeV and other similar laboratory tests do not 
constrain the OP model. Provided < 1.3 x 10~^^eV 
in the galaxy, and it was shown in ijllBI that this is 
expected to be the case, the starlight polarization con- 
straint on the OP model translates to: 

A., 1/2 



G^'^Mq > 1.6 X lO^TeV 



10- 



(53) 



at 99.9% confidence where da/ a is the fractional differ- 
ence between a in the laboratory and a in a background 
such as the galaxy. For comparison, the previous best 

constraints were Afo > 15TeV and ^p^^'^Mo > 3TeV. 
If |(5a/a| ~ 0(10"^) as suggested by the analysis of 
Webb et al. [29| . then this represents an improve- 
ment of two to three orders of magnitude. We note 
that if a subsequent analysis were to confirm the lower 
bound on 1/Af found from starlight polarization mea- 
surements, then both these measurements and the Webb 
et al. value of 5a/ a could be explained by an OP model 
with 



-1/2 



Afn 



8 X lO-' TeV. 



and Al - (0(10-2) - 0(10^)) eV; 0„ - (10 - 20) TeV. 
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VII. CIRCULAR POLARIZATION: A SMOKING 
GUN? 

The total polarization due to chameleon-photon mix- 
ing grows as the square of the frequency of the light 
until it reaches a critical frequency at which the mix- 
ing becomes maximal. It has not been possible to de- 
tect this frequency pattern in current linear polariza- 
tion data. Objects whose initial polarization is well con- 
strained have not been observed over a wide enough fre- 
quency range or to the required accuracy to see any such 
signal. Certain GRBs have been observed over a very 
large range of frequencies as they evolve, but because 
their initial linear polarization is not known accurately, 
and generally does not satisfy po Ij it is difficult to 
search for a chameleon signal in this data. 

The production of circular polarization by chameleon 
photon mixing has a much more interesting signature. 
One does not usually expect significant amounts of intrin- 
sic circular polarization (CP) for astrophysical objects. 
We noted above in mV CI that chameleonic CP produc- 
tion is peaked over a potentially large range of wave- 
lengths (when iV > 1) i.e. Aosc = Krit/N < A < A,„ax, 
where 



Simulated Data for Light Passing through a Gaiaxy Ciuster 



Acrit max 1, 



nbl\ 



ttM I 



Importantly, in this band, the chameleon induced cir- 
cular polarization is the same order of magnitude as the 
chameleon produced linear polarization, and both exhibit 
a highly oscillatory frequency dependence in this region. 
Outside of this wavelength band, the chameleon contri- 
bution to the circular polarization is much smaller than 
to the linear polarization. If mixing is maximal, q ~ 0(1) 
is expected. Neither the magnitude, the shape, nor the 
oscillatory frequency dependence of the chameleon in- 
duced circular polarization peak is likely to caused by 
any other process. The observation of this peak could be 
considered a smoking gun for chameleon-photon mixing, 
and if such a structure could be ruled out then strong 
constraints on chameleon like theories would follow. In 
particular if 0(1), highly frequency dependent values of 
q are not seen in the region Aosc — Acrit /.^ < A < Acrit, 
maximal mixing could be ruled out, immediately limiting 
M > 10^° ~ 10" GeV. Strong constraints would result if 
the CP of a distant object whose light was known to pass 
through the magnetic field of a galaxy cluster could be 
constrained in the region Aosc < A < Acrit . To ensure the 
maximal sensitivity to chameleonic effects however the 
spectral resolution would have to be w Aosc or smaller, 
which for a cluster would require SX < O.lA. Assum- 
ing light travels roughly 1 kpc through the galaxy, 1 Mpc 
through a galaxy cluster and about 2.5 Gpc through the 
IGM the typical expected values of Aosc and Acrit are 
shown below in Table |lTl We have assumed ^ 
6.4 X 10-12 eV in the galaxy, <C 2.2 x 10'^^ eV in the 
ICM and < 1.8 x IQ-i^eV in the IGM. 




FIG. 5: Simulated data for two objects whose light has passed 
through roughly 1 Mpc of the magnetic field of a typical 
galaxy cluster. We have assumed <g 2.2 x 10"^^ eV and 
M — GeV; which corresponds to strong mixing for wave- 
lengths of Amax ~ 24Acrit . We have assumed that both objects 
have little or no intrinsic circular polarization, and are 50% 
linearly polarized prior to chameleon mixing. Qualitatively 
similar behaviour is seen for different values of the intrinsic 
linear polarization, mjo, and in particular the behaviour CP 
fraction does not depend greatly on mio- 



FIG. shows simulated data for two objects (e.g. 
GRBs), with 50% initial linear polarization and no intrin- 
sic circular polarization, whose light has passed through 
about IMpc of the magnetic field of a galaxy cluster. 
The wavelength. A, in this plot should be interpreted as 
'^m/ (1 + ^ciust) where Am is the measured wavelength and 
^^ciust is the redshift of the cluster. We have assumed, as 
is generally the case, that ^ 2.2 x 10--^^ eV. We 
also have taken M — 10^° GeV, which corresponds to 
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TABLE II: Position of CP Peak 



Environment 




Galaxy 


I2A 608A 


ICM 


0.1 A 104A 


IGM 


1.5 A 3600 A 



Simulated Data for Light Passing ttirougti 1 kpc ot Itie Galaxy 




in the region Aosc < A < Acrit, potentially detectable lev- 
els of CP (between 0.1% and 0.5%) are typical. 

If < 9 X 10— 12eV, measurements of CP between 
A 0(1)A and A ^ O(1000)A for astrophysical objects 
should allow one to detect or rule out theories with M < 
lOiOQeV. 

Thus far, circular polarization has been measured for 
a number of different astronomical sources; for certain 
stars observed in the near infrared in [s^, [sj] , for zodia- 
cal light in jl^ , for the Orion molecular cloud in ^sd] , for 
some relativistic jet sources at radio wavelengths in ST}. 
However for all of these observations mixing between pho- 
tons and chameleons is weak, and the frequency resolu- 
tion of the observations is not good enough to detect 
a chameleon signal. Additionally all such observations 
have been made at wavelengths outside the expected 
Aosc — Amax positlon of any chameleonic CP peak. 



VIII. SUMMARY 



FIG. 6: Simulated data for two objects whose light has passed 
through roughly 1 kpc of our galaxies magnetic field. We have 
assumed < 6.4 x 10"^^ eV and M = 10^° GeV. We have 
assumed that both objects have little or no intrinsic circu- 
lar polarization. Potentially detectable levels of CP are seen 
between Aosc ~ 12Aand Acrit 



strong mixing for A < Amaxi in this case Acrit ~ 104A. 
We can see that chameleonic production of polarization 
begins for A < Amax and here Amax ~ 24Acrit ~ 2500A 
(i.e. in the middle UV part of the spectrum). Very simi- 
lar behaviour is seen for different choices of the intrinsic 
polarization. Between Aosc and Amax both the linear po- 
larization, mi, and the circular polarization, rric, are, as 
expected, highly frequency dependent and as we expect 
from the strong mixing scenario when Aosc < A < Amax 
the magnitude of both rric and m; oscillates between 0% 
and 100%. For A < Aosc, mi/100% settles to some, essen- 
tially random value between 0% and 100%, and nic 0. 
If there is little or no intrinsic circular polarization, the 
behaviour of tUc does not depend greatly on the value of 
the intrinsic linear polarization. If such measurements 
could be made it should be straightforward to either 
detect or rule out values of M for which strong mix- 
ing in clusters is expected to occur; M < 10^^ GeV. 
All light that reaches us from distant objects will have 
passed through at least a part (~ O(l)kpc) of our own 
galaxy's magnetic field. In FIG. [6l we show sample cir- 
cular polarization data for two objects (with little or 
no intrinsic circular polarization) whose hght has passed 
through 1 kpc of our own galaxy's magnetic field (corre- 
sponding to about 50 regions of the random field, and 
one of the regular field) . Most of the CP production is 
due to the random field. We take M = 10^° GeV and 
<C 6.4 X 10~^^eV. In this case we are in the weak 
mixing limit, and Amax — ^mt = 608 A. We can see that 



Theories of physics beyond the standard model typi- 
cally predict the existence of new scalar fields. If these 
scalar fields do exist it is important to understand both 
their self interactions and their interactions with the 
other fields present in the model in order to test and 
constrain the theory. In this article we have studied the 
results of a coupling between scalar fields and photons 
on observations of astrophysical objects. Specifically we 
have studied the scalar fields of the chameleon and Olive- 
Pospelov models, which are strongly interacting in low 
density environments yet currently undetected in the lab- 
oratory. For simplicity we refer to both types of scalar 
field as chameleons. 

If the chameleon field couples to photons then in the 
presence of a background magnetic field the chameleon 
mixes with the component of the photon polarized or- 
thogonally to the direction of the magnetic field. We have 
studied the effect of this mixing on light beams passing 
through a large number of randomly oriented homoge- 
neous magnetic domains, in order to predict the effects 
of chameleon-photon mixing on observations of light from 
astrophysical objects. Typically both linear and circular 
polarization are induced in the light beam by mixing in 
such an environment. 

We found analytic solutions to the equations describ- 
ing the mixing in two important limits. In the weak mix- 
ing limit the polarization fractions induced by chameleon 
photon mixing are highly wavelength dependent. If the 
light is not polarized at the source the averaged values of 
the total and circular polarization scale as NP^^^, that 
is as the product of the number of domains traversed 
and the probability of mixing in any one domain. This 
limit is generally appropriate when one is considering the 
chameleon induced polarizations at wavelengths longer 
than roughly lOOOA. In the maximal mixing limit, which 
applies when the chameleon-photon coupling is strong 
and the wavelength is sufficiently short, little or no cir- 
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cular polarization is produced by the mixing, but the 
production of Hnear polarization is at its strongest. The 
distribution of the total polarization fraction after mix- 
ing in a large number of domains is independent of the 
parameters of the chameleon model, and instead depends 
only on the initial polarization of the light. The average 
value of the total polarization fraction is always greater 
than (7r/2) — 1 w 0.57 in the maximal mixing limit. 

Numerical simulations confirm the analytic analysis. 
In particular they clearly demonstrate the existence of 
two wavelength scales Aosc = Krit/N = 47r2/|m^g|Lpath, 

and Aniax = Adit max (^^, tt^ai ) ^^^'^^ determine the 
shape of the polarization signal. Here B and L are 
the strength and domain size of the magnetic field. 
3077 = is the coupling between two photons and 

the scalar field. The linear polarization is greatest for 
A < Amax, and the circular polarization is peaked for 
Aosc ^ A < Amax. Both polarizatiou fractions are highly 
frequency dependent for A > Aosc- This highly oscilla- 
tory behaviour means that observations of polarization at 
these wavelengths must be performed with a sufficiently 
good spectral resolution if any chameleon induced signal 
is to be resolved. 

We have considered a wide variety of astrophysical 
observations and have used these to constrain the pa- 
rameters of the chameleon model. From observations 
of starlight polarization in our galaxy we show that at 
the 99% confidence level M > 1.1 x 10^ GeV, which 
is an improvement of over two orders of magnitude on 
the previous best constraints on M . The equivalent con- 
straint on the Olive- Pospelov model is given in ([55]) . Both 
constraints could, however, be evaded if the potential of 
the scalar field, V{(f>), is chosen so that the field is suffi- 
ciently heavy in regions with the density of our galaxy: 
10~^^ eV. Constraints from objects outside the 
galaxy are limited by our lack of knowledge about a 
possible intergalactic magnetic field, -Bigm- If, however, 
-SiGM ~ 10~^ G and is coherent over roughly Mpc scales, 
then the lower bounds on M and the OP model coupling 
scale are raised by roughly two orders of magnitude. 

The circular polarization signal predicted from 
chameleon-photon mixing in a large number of randomly 
oriented magnetic domains was shown to be very distinc- 
tive. Its frequency dependence is unlikely to have been 
caused by any other physical process, particularly as as- 
trophysical objects do not normally produce significant 
amounts of circular polarization. To date, no observa- 
tions of astrophysical circular polarization have yet been 
made with sufficient accuracy to allow us to search for 
a chameleon signal. Nonetheless, we have shown how 
future observations of circular polarization in the wave- 
length range 0(1) — O(1000)A would be a smoking gun 
for chameleon-photon coupling. 

We have also reported a seemingly strong statistical 
preference in observations of starlight polarization in our 
galaxy for the presence of a chameleon-like field. Pre- 



cisely, at the 99% confidence level, we find 

where B and L are the strength and domain size of the 
random component of the galactic magnetic field. For- 
mally, the central value deviates from zero (the value for 
a theory without a chameleon) by more than IOct. It 
must be stressed, however, that this is only a prelimi- 
nary analysis and we have only performed it for three 
out of a possible 121 objects. Before a detection could be 
claimed with any confidence, a full study of the possible 
backgrounds and systematics for these observations that 
could bias one towards larger values of 1/M would have 
to be undertaken. Based on initial numerical simulations 
of data, it does, however, seem unlikely that polariza- 
tion due to interstellar dust would produce such a strong 
signal. 

In summary: astrophysical polarization measurements 
currently provide the strongest constraints on any cou- 
pling between photons and the scalar field, for many 
chameleon and chameleon-like theories such as the Olive- 
Pospelov model, improving on previous constraints by 
more than two orders of magnitude. Furthermore, fu- 
ture measurements of linear and, in particular, circular 
polarization at short wavelengths (i.e. < 2000A) could 
provide one of the best tools in the continuing search for 
such scalar fields. 
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APPENDIX A: FLUCTUATING ELECTRON 
DENSITY 

In Ref. [131 it was shown that fluctuations in the elec- 
tron density, Uo, and hence the plasma frequency, Wpi, 
could lead to a significant enhancement of the photon to 
axion-like-particle (ALP) conversion rate when ^ ujpi 
and |A| « LUp^L/Auj ^ 1. In this appendix we reproduce 
this analysis and show that conclusions about the mag- 
nitude of such an enhancement effect are modified in the 
light of more recent models of the electron-density in our 
galaxy (specifically the NE2001 model) than those used 
in Ref. [4j. We also extend the analysis to allow for 
fluctuations in the magnetic field B. 

We write a;pj(z) — + Sn{z)) where Sn{z) = 

Sne{z)/ne. We also have B = \B±\ = B{1 + db{z)). 
We are concerned with the limit to^ <C Wpj. Now in a 
single region of magnetic field the equations describing 
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the evolution of the photon and the chameleon are: 

- 711 +711,2^ = ^l\{zh\\, 



B 



B 

-4> + (l),ZZ = --^7-1.2- 

We assume that uj ^ ujp\ and 5'„/6n ^ us- We write the 
solution for 7|| thus: 

7|| = 7oe ^ ' ^ ' ■ 

where we assume Wpi <C oj and so \a_z \ <C lo. We similarly 
write: 

0(^) = 0(z)e"^(^-*)-*°("). 



We then have: 



^0, 



B 



lUJ, 



We define m^g = — cZ'pi(l + ~ — Wpi(l + 5„) where 
z5n{z) = / (5„(z') dz'. 



We then let a; = m^fjz/4a;, so z = L implies x = A and 
remember that: 



where 



tan 26* 



'UjI^. Thus: 



u 



- 7 



2Boj 



iM(x) 



where 



M{x) 



I -iD{x)ia,n2e 

y iD(a;)tan26l -2 

(ct3 - I) + fT2i:'(a;) tan 20. 



(Al) 



(A2) 



where 



D{x) = l + 5{x) = ^ ' 



In the above the ai are the Pauli matrices. We write 
M{x) = Mo + cr2^(a:)tan26' where Mq = (^3 - I) + 
CT2 tan26' and so: 



.sf^) 
V cos 261/ 

+j((73 cos 29 + (72 sin 26) sin 



cos 26*. 



where 



^(x) = sin 20 sin 
tan<p(x) = cos 20 tan 



cos 20/ ' 

X 



.cos 20/ 
We note that 

g-^Mox^^^^Mox ^ (^os 20ni + sin20n2) • a (A3) 

where 



ni(x) 



;os20cos (-^) 

V cos 29 / 

sin 20 cos I 




cos 2e ; 



n2(a;) = | sin 20 
cos 20 

We then define v = g-«c(x)g-iMo2;^ -^y^gj-e 

C(x) — (n2 • a) tan 20 / d{x )dx = — . 

/n cos 20 



We then have: 

v^x = i sin 295{x) (n{x) ■ <t) v. 

where 

COS 20 COS ^ „; 

y cos2^? J 



(A4) 



n 



V 



sin 20 cos 



/ 2(l+a)a: \ 
\^ cos 29 J 



I 



When it is acceptable to do so we may solve Eq. (|A4p 
perturbatively, a sufficient condition is: 

||Asin20(5(a;)|| < f, 

as r]?{x) = 1. This may be satisfied if either 20 <C f, 
A <C f or ||(5|| <C 1. To sub-leading order we have: 

V w N(x)vo = [I + isin20(c(j;;0) -cr) 



(A5) 



where 



isin^ 20 (c^ + id-a) 



c{x-e) = I ds(5(s)n(s;0), 



vo, 



d{x; 



As~5{s) (n(s;0) x c(s;0)). 



(A6) 
(A7) 



We evaluate N(x) in the weak-mixing limit, 20 ^ I. To 
do this we expand the diagonal terms in N(x) to order 
(20)2 and the off-diagonal ones to order 20. We find 



N: 



I-202(||.g||2+ir) 
-205 



20g* 

I-202(||g||2-^r) 



,(A8) 
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where: 



Thus: 



T = 2Re(g) + h, 



h = d3(x). 

Now at the end of a magnetic domain with length L, 
z — L and u — u{L) we have: 



u(L) 



-iA 



(l + a(£))A \ 

cos 20 y 



cos 

(l + a(L))A\ 



(A9) 



cos 26* 

(era cos 29 + (72 sin 26*)] N(i)u(0) 

where u(0) is the initial value of u. In the weak mixing- 
limit, it follows that the probability of converting a pho- 
ton to a chameleon is: 



46*2 1 1 sin 



(l + a(A))A\ 



i29 



1 -5(x = A)||^(AlO) 



The second term inside the || • |p represents the en- 
hancement term from electron density fluctuations. It 
is straightforward to check that if 20 <C 1 we expect 
g/A < 0(1) if, as expected, H^f < 0(1). Therefore 
when A ^ 1, i.e. at high frequencies, we do not ex- 
pect the new term to produce a large enhancement. This 
was also noted by Carlson and Garretson in Ref. [i^ . 
We therefore focus on the limit A ^ 1. In this limit 
sin^((l -I- Q!)A/cos26') ^ 1/2 on average. Wc denote 
the relative magnitude of the enhancement in photon to 
chameleon conversion by E and 

E = 2\\gr. 

When E ^ 1, the enhancement is negligible, and when 
-E 3> 1 the enhancement is strong and represents a signif- 
icant effect. In the weak-mixing limit in which we work 
llall ^ 1, and so: 



E ^ 2\\ I da;J(x)e"f 
/o 

Jo 



(All) 



We define the Fourier transform Sk(k) of S{x) thus 
S{x) =5{Az/L) ^ I d3fc4(k)e''=^^^/'^. 



where fc^ is the component of k in the z direction. The 
power spectrum P{k) is given the by expectation of 
4(k)J^.(q) thus: 



4(k)^:(q)\=P(fc)<5(3)(k-q). 



{E) = 2A^ J d^kP{k)smc^ 



k^L + A 



(A12) 



where sinc(a;) = sinx/x. Electron-density and magni- 
tude field fluctuations are often modeled by a Power spec- 
trum with inner scale Iq and outer scale Lq, and a power 
law behaviour between these two scales i.e.: 



P(fc) 



C2 



1 a/2 



(A13) 



When, as is the case for visible light, AIq/Lq <C 1, the 
role of /q in the estimate for E is negligible, and we may 
approximate by setting la — Q- We then find that: 



E=2> 



where A^ 



ds5{s)e' 



S^^A^C^Ag-" 
(a - 2)L 



, (A14) 



A^L ^. It can be similarly checked 



that with this form of P(k): 



and so 



E^iqA^\r5' 



1 + A2L2/L2 



(A15) 



(A16) 



where 



3(a-2)r(f-§) 



For a = 11/3, which corresponds to a Kolmogorov power 
spectrum, qa ~ 0.9958.^ 

We remember that 5 — {5b{z) — 5n{z))/{l + (5„(z)) 
where ^j, is the magnetic field fluctuation and (5„ is the 
electron density fluctuation. The power spectrums of 
both fluctuations are generally taken to be described by a 
Kolomogorov power spectrum with some inner and outer 
scale. We note that if, as is often assumed, the inner 
and outer scales of the magnetic and electron density 
fluctuations are the same, and if the two fluctuations are 
uncorrelated when i5„ <C 1, (5 will also have a Kolmogorov 
type power spectrum. We also note that correlations be- 
tween 5b and (5„ could potentially greatly decrease the 
power in 5. Specifically if Si, w (5„ then 5 Sb,Sn. The 
structure of electron density fiuctuations in our galaxy is 
much better understood than the structure of magnetic 
field fiuctuations. For simplicity, and to make an order 
of magnitude estimate of E we take Sb — and assume 
6n <C 1 so that S ~ —Sn- The power spectra of (5„ and S 
are then equivalent. For electron density fiuctuations, es- 
timates of the inner scale, Iq, place it around 10^—10^, m. 
It can be checked that for ui > 10~^ eV and L w 50, pc, 
AIq/L <C 1 as assumed above. 
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In the NE2001 model [48[ for galactic electron density 
fluctuations, the fluctuation parameter 
thus: 



Fn, is defined 



1 pc 



2/3 



and the electron density fluctuations have a Kolmogorov 
spectrum with a ~ 11/3. We also estimated previously 
that: 



16 



2eV 



and so 



E « 15AF„ 



La 
1 pc 



where 



/3 = LlA~- 



1 + 0.1 



5/3 



1 pc 



-5/6 



(A17) 



2eV 



The fluctuation parameter varies widely across the 
galaxy. On average in the disk ~ 0.2 however in 
the local interstellar medium (out to about a kpc from 
the Sun) F„ « 0.01 — 0.1. The stellar objects we ana- 
lyzed in WI Al are located in the local ISM where F„ is 
smaller, however even if we take the slightly larger value 
of Fn ~ 0.2 appropriate for the disk on average, we find 
that for visible light w ~ 2 eV: 



E ■ 



Lq 
1 pc 



5/3 



-5/6 



Carlson and Garretson [44| took the outer scale of tur- 
bulence to be Lq « 10 — 100 pc, which results in E be- 
coming independent o f Ln and E ^ 19 — 20. However, 
more recent estimates [89| suggest a much smaller value 
for Lq than previously expected, specifically an Lq that 
is no more than a few parsecs. In HII regions (clouds of 
gas and plasma in which star formation is takii ig p lace) 
Lq « 0.01 pc and the pulsar measurements [90| give 
Lq w 0.03 pc. The precise value of E therefore depends 
fairly strongly on the value of outer scale for the Kolo- 
mogorov spectrum Lq, which is uncertain. This is be- 
cause the enhancement term is predominately sourced by 
electron fluctuations on scales of ^ ~ L/A « 4ti;/||m^jf || . 
For visible light, I ~ 0(l)pc. The structure of galactic 
electron density fluctuations is not, however, well under- 
stood on such scales and almost all measurements of such 
fluctuations relate to lower scales. This means it is dif- 
flcult to make an accurate estimate of the enhancement 
factor. However given Lq < fewpc, lu ~ 2,eV, we esti- 
mate 0.03 ^ / < 10 based on the different estimates for 
Lq. a correlation between electron and magnetic fluc- 
tuations could signiflcantly lower this estimate. Hence 
we have estimated i? to be 0(1) or smaller in the visible 
part of the electromagnetic spectrum. Importantly, even 



if the conversion rate is enhanced, the oscillatory nature 
of the chameleon induced polarization remains. Thus a 
slightly enhanced conversion probability is not expected 
to signiflcantly alter the form of the signal for which we 
have searched. Given the great ambiguity in the precise 
magnitude of the enhancement term and because we es- 
timate it to be no greater than factor of about 10, we 
have chosen to neglect it in our analysis. 

We now consider the magnitude of any enhancement 
effect due to electron density fluctuations in galaxy super- 
clusters, such as that considered by Jain et al. in Rcf. 
[4^ . Since very little is known about electron density 
fluctuations in galaxy clusters and superclusters, Jain et 
al. assumed a simple scaling relation where all unknown 
dimensionful quantities scale with rig. They did not how- 
ever include the role of an outer scale of fluctuations, Lq, 
instead assuming that P{k) was everywhere a power law. 
The outer scale of fluctuations is important as it is re- 
quired for the total magnitude of fluctuations (5^) to 
be flnite. We assume the same scaling for dimensionful 
quantities as that used by Jain et al. . We therefore 
assume that the length scale Lq cx n7'^^^, but that ^(S^) 
is approximately the same in a galaxy cluster as it is in 
the galaxy. 

If Lq — Ipc in the galaxy where rig ~ 0.03 cm~^ then in 
the galaxy supercluster considered by Jain et al. where 
He « 10~^cm~'^, one would expect Lq ~ 31 pc. In the 
same region tjpj « 3.7 x 10^^^ eV and an appropriate 
value for L, the length of the magnitude domain, is sug- 
gested in Ref. [i^ to be 100 kpc. This gives: 



A w 2.7 



2eV 



(A18) 



Thus using Eq. lAlGl for a = 11/3 and assuming i^S^J^ 
1, we flnd that for a galaxy supercluster the enhancement 
factor for visible light is estimated to be: 



E < O{10-^) 



31 pc 



and so any enhancement due to electron density fluctua- 
tions in this region is estimated to be sub-leading order. 
Jain et al. found the opposite result but ignored the role 
of the outer scale, Lq, which limits the overall magnitude 
of fluctuations. 



APPENDIX B: CHAMELEON OPTICS FOR 
MULTIPLE MAGNETIC DOMAINS 

In many realistic astrophysical settings, light beams 
pass through many magnetized domains, and in each 
domain the angle of the magnetic fleld relative to the 
direction of propagation is essentially random. In this 
appendix we present, in detail, the equations which de- 
scribe this multiple domain problem and their solutions 
in a number of important limits. In ijlV Al we presented 
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the equations that describe how the chameleon and pho- 
ton fields evolve as they pass through a single magnetic 
domain. In that section we split the photon field into 
components polarized parallel and perpendicular to the 
direction of the magnetic field, and used this as a basis 
to define the Stokes vector, (/-y, Q, U, V)'^ for the photon 
field as well as four associated amplitudes, J, K, L and 
M, which describe correlations between the chameleon 
field and components of the photon fields (see Eq. (fT6|) 
for the definition of these quantities). To deal with the 
multiple domain case we must first fix a basis for the pho- 
ton field that is independent of the direction of B. Doing 
this we take the two components of the photon field to 
be 71 and 72, and redefine 

A = (l7iP) + (l72p), 

Q = (i72n-(i7ip), 

2 (7271) , 
2e^'^ (7iX> , 
^e'" (72X) , 



U + iV 
J + iK 
L + iM 



and as in ijlV Al we define X = 3/^ — 2. 
independent of the choice of basis, 
in the n*'' magnetic domain: 



and V are 
We define 6n so that 



and define 



7|| = cos6'„7i - sin6'„72, 
7j_ = cos6'„72 -f sin6i„7i. 



[/' 
J' 

L' 
K' 
M' 



(3cos26'„ + ;7sin26'„, 
-(5sin26'„ + C/cos26'„, 
J cos 9n — L sin 0„ , 
J sin 9n + L cos 6n , 
K cos On — M sin 6'„ , 
K sin On + M cos On ■ 



V in the -"^ 



The evolution of the primed quantities as well as X and 
^th j-ggjon are then described by Eqs. (fT?]) - 
with Q being replaced by Q', U by U' and so on. 
Solving the full system of equations for iV 3> 1 do- 
mains involves diagonalising an 8 by 8 matrix as well 
as evaluating multiple sums involving the random angles 
0„ for n = to — 1, and we have been unable to 
find an analytic general solution. It is straightforward 
to solve the system numerically, but analytical solutions 
are often more useful for understanding the behaviour. 
Fortunately, it is possible to make a great deal of ana- 
lytical progress in the weak-mixing limit where Na <ti 1 
and NPj^^ <^ 1, where N is the number of magnetic 
domains, as well as in the strong mixing limit where 
NA < 1 and NP^^^ > 1. 



these limits ip Ki A, (3 ~ 2 A and 
tan^ 261 



a = Lp — A: 



2A- sin2Al 



(Bl) 



We assume that there is no initial chameleon flux, so that 
initially X ^ Xq = I and J = K = L = M = 0. Without 
loss of generality we pick our coordinate basis so that 
Q — initially and U = Uq and V = Vq. By requiring 
that NPj^^ <^ 1 and Na <C 1, we are assuming the 
perturbations, SX, SQ, 6U and 6V, are small compared 
to the quantities X, Q, U and V, and that J, K, L and M 
are <^ 1. We define J — Aj and make similar definitions 
for fc, m and n. We compute the perturbed quantities to 
OiNA^) and 0{Na^). We define (5X„ to be the value of 
SX after having passed through the n*'' region, and make 
similar definitions for the other quantities. Expanding 
to first order in the perturbations we find the following 
simplified recurrence relations 



Q /12 oa2 _ 

SXn+l = SXn-— —UQSin20n (B2) 

+3A'^ {In COS On + jn sm On) sin 2 A 
—3A^ {nin COS On + kn sin On ) COS 2 A , , 

SQn+1 - SQn~^COs20n (B3) 
+A^ {I n COS Ufi Jn sm C7yj ) sin 2A 



n COS Ufi Jn sm Un ) sm ZZA 
—A^ {rUn cos On — kn sin On) COS 2 A, 



—aVn sin 20n + — sin AOnUi 



SUn+l = SU, 



A- 



n 

VA^ {In sin On + in COS On) siu 2 A 
-A^ {nin sin On + kn COS On) cos 2A 

Q,2 

^aVn cos 20n —Uq COS^ 20n 



A^ 

— sm26'„ - -^C^o 



(B4) 



Wn+X = SVn - ^Vo 

+A^ {In sin On ~ jn COS On) COS 2 A 

+A^ {nin sin On — kn COS On) sin 2 A 

1 2 

--a Vb - aUn COS 26'„, 
+aQn sin 20n 



(B5) 



and 



1. Weak Mixing Limit 



When iVa < 1 and NP. 



7^0 



we 



must have either 

A/ COS 26*, A tan 2 A < 1 or tan 261, Atan^26l < 1. In 



kn+l + ijn+1 

rrin+i + iln+1 

Yn 



.2.A, 
„2iA , 



6" {kn + ijn) + Yn, 



(B6) 

e— (m„ + iln) + Zn, (B7) 
{Uo + iVo) cos On + sin On, 
{Uq — iVo) sin On + cos On- 
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Eqs. (|B6D and ([BTD are solved thus 

n-l 

kn + i]n = ^ (B8) 

n-l 



and 



r=0 



Assuming 3> 1, we then arrive at following solutions 
for the perturbations to the components of the Stokes 
vector to 0{NA^, Na^): 



SXn = -■ 



3NA^N 



(B9) 

6Qn = -AA2i?^+'(2A)- AA2[/oi?^"^(2A) (BIO) 
(2iVA2 + Ara2^ 



SVi 



N 



4 -C/o (Bll) 

+NA^VoQ^^'^\2A) + VNaVon'k 

(NA^ + Na^) f^_^, , 
^ ^-Vo-NA^e^^ '(2A) (B12) 

-AA2c/ogi^+'(2A) 
-A^A2yoi?^"^(2A) - VNaUoK% 



where 



AT-l n-l 



1 " 

^?^±(2A) = - ^ ^cos(2A(n-r))cos(0,±0„), 

n=0 r=0 
Af-ln-1 

(2A) = ^ ^ 51 cos(2A(n - r)) sin(0, ± 0„), 

n=0 r=0 
Af-ln-1 

e^±(2A) - ^ 5] ^ sin(2A(,i - r)) cos(a, ± 0„), 

n=0 r=0 
AT-ln-l 

g'j^ (2A) = ^ 51 E sin(2A(n - r)) sin(0, ± 0„). 



JV-l n-l 



m5v 



cos 26'„ cos 2^r 



TV E E' 

n=0 r=0 
Af-ln-1 

COS 20 J., 

n=0 r=0 
Af-ln-1 

— 5] sin26'n sin 26>r, 



N 

n=0 r=0 

-^5:cos20„, 

^ n— 
1 JV-l 

n=0 



(B13) 



Each of these nine quantities vanishes when averaged over 
all possible values of When f{n,p) = f{p, n) we have 



N 



n=0 p=0 



n=0 p=0 
N-l 



-2^E/("''^)' 



and so 



,qc± _ 



n=0 



\ (xl T Xl + Xl T - i ± i ) (B15) 



where 



XccXcs + AscAss, 
X scX cs Xg^X^Q. 



(B16) 
(B17) 



JV-l 



A, 



Xr 



X, 



X.. 



— cos 2nA cos 6'„, 

n— 
1 JV-l 

^ > cos2nAsin0„, 

n— 
1 JV-l 

— ^ >^ sin 2n A cos dn 

^ n— 
1 JV-l 

> sin 2n A sin 0„ . 



n=0 



In the large N limit (and at fixed A) each of these four 
quantities are independent, normally distributed ran- 
dom variables: X^c, X^s ~ ^(0,(7^) and X^^, X^s ~ 
7V(0,cr^), where: 



4 



1 ± 



cos(2(iV- l)A)sin2iVA 



iVsin2A 



(B18) 



Additionally, i^^, fp^, are, for fixed A and in the 
large A^ limit, well approximated by independent nor- 
mally distributed random variables, with ~ A^(0, <j\) 
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and the rest are n(0, cr|) where for iV ^ 1 



sin^ 2A^A 
iV2lh?2A 

sin^ 2A^A 
iV2 sm^ 2A 



We choose a basis so that initiaUy L/q > and define 
m;o = f/o, TOco = Vb, go = l"icol and po = VC/q + Vq^. 
Keeping terms to order 0{N P^^fppo) and 0{N^P^^^), 
we find that: 

+N^P^^^{l-plf(^l+r^^ . (B19) 

where we have used Eqs. (jB15IIB17| to provide the iden- 
tity: 

1 



,(c+)2 

'n 



'n 



{s-)2 



Xi 



N 



(B20) 



If Po ^ 0(1), the last term in this expression is the same 
order as terms that have been omitted so it too should be 
dropped. Similarly for the fractional circular polarization 
we have to 0{NP^^^) and OiNa^): 



mc{N) 



mcQ 



(B21) 



-iVP^^^(l-m2o)e^-(2A) 

+NP^^^miomcoi}'+{2A) 

-NP^^^miag'^{2A). 

If there is no initial polarization (po = 0), or more 
generally if NPj^if,(l — p§)/po S> 1, the final polarization 
fraction is given by 



p{N) 



NP^ 



(B22) 



We may therefore write 
1 



p{N) 



-NP 



aliXl + Xi) + a^_iXi + Xj)) , 



where the Xi are independent identically distributed 
A''(0, 1) random variables. When po = the circular 
polarization simplifies: 

mc{N) = NP^^^a+a^ {XiX^ - X2X4) . 

Where po 7^ and NPj^^{l — Po)/po <C 1 we have to 
0(7VP-,„^(l-p2)/po): 



p{N) 



Po 



NP^^^{1 -pl)mio , 2 



-a'^_X'^Xj^ 
NP^^^{l-p'^)mco 
Po 



a,XiX2 



a+a- {X1X3 — X2X4 



The circular polarization is given by Eq. (|B2T1) in this 
case. 



2. Strong Mixing Limit 

We now consider the strong mixing limit. This is the 
limit in which iVA ^ 1 so that P^^^f, takes it largest 
value, and the mixing between the chameleon and pho- 
tons is strong, NP^^^ 3> 1. In this limit a, /?, A, (/3 <C 1 
and so Eq. (|17ll22p simplify to: 



X 

Q 

U 
M 
K 



1 - -A^\ X - -AH 
2 2 



~3AVi -A'^M, 



1 



(B23) 



-A\/l - A^M, , 

Vl - A^U - AK, 

(1 - 2A^)M + A^1~A^Q + X). 

Vl - A^K + AU, 



and 



V 
J 
L 



^/l-A^V - AJ, 
\/l - A'^J + AV, 
L. 



The diff'erently oriented magnetic fields in each domain 
mix Q with U, M with K and J with L. It is clear then 
that the evolution of J and L are completely decou- 
pled from that of X, Q, U, M and K. We are concerned 
with the limiting value of total polarization fraction, p. 
Additionally since we expect the initial circular polariza- 
tion fraction to be small, qQ = |mco| ^ Po, we set V — 0. 
We also require that initially the chameleon flux is zero 
(M = K = L= J — initially) . It is clear then from the 
above equations that V remains zero. From simulations 
we see that in the strong mixing limit the final mean po- 
larization fraction takes a specific value, which depends 
on pq. Remarkably we can calculate both the limiting 
value and the final distribution of p analytically without 
actually explicitly solving the above equations. 

We assume that initially the photon is in a state with 
polarization fraction po = (1 — a)/(l + a). Without loss 
of generality we pick coordinates so that U — initially 
and write the initial Stokes vector of the photon state 
thus: 



Sq 



('A 






Q 




(1-a) 


u 







\v) 




I / 



(B24) 



We can always consider such a partially polarized photon 
state to be a linear superposition of two fully polarized 
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photon states (labelled (+) and (— )), i.e. 

Sa = S+{Q) + aS-{0) 
where (dropping the V component as it vanishes): 



S±{0) 




(B25) 



Since both S'+(0) and 5'-(0) represent fully polarized 
photon states, they can also be described in terms of a 
vector whose components are the photon and chameleon 
amplitudes, Ci = 71, C2 = 72 and = x = We define 
this vector to be w+ for 5+ and u_ for , so that 



ci = 1 • (B26) 




(B27) 



We also define Wtot = w+ + w- and note that this too is a 
fully polarized state. 

The evolution of a fully polarized state through a single 
magnetic domain is given by Eqs. (|12m3p . We note that 
these equations conserve the total flux 1^ + ~ Aq, 
where = |cip + |c2p and = jc^p. For v±, the total 
flux is 1 and for Utot it is 2. 

After having passed through many randomly orien- 
tated magnetic domains, if NPy^^ :§> 1, the mixing 
between the chameleon and photon fields, and between 
different components of the photon field, will be strong. 
This means that on average the initial flux should be 
evenly distributed among each of ci, C2 and and so 







1 


Cl 






V^2. 


N 





\J\ — cos Q 



where each of c^, ci, C2 are uniformly distributed random 
variables on Ao[— 1,1). This implies that x ~ t7[— 1,1) 
and Q - [/[0,27r). 

Now u_ and Vtot are all fully polarized states. If, 
after having passed through many regions, w+ w+(oo) 
and u_ V-{oo) where 



«+(oo) = 



u_ (00) 



y/1 — COS 9 
y \/l " sin 9 



cos ( 



^l-?;2 sm( 



then since the field equations are linear i>tot — > vtot{oo) = 
w+(cxd) + w_(cxd). Now in the limit of strong mixing, the 
C0 components of u+(oo), w_(oo) and vtot{oo) must all 
be uniformly distributed random variables on [—Aq,Aq). 
This imposes a very strong condition on the distributions 
of X and y, in fact one must have x Vl — X"^ cos ■0 and 
y = Vl — sinip where ip and X are independent uni- 
form random variables: ^ ^ U[0, 2tt) and X ~ U[0, 1). 
We also know the total flux, utot- Initially the to- 
2, 



tal flux is Aq = 



and finally it is Aj — {x 



y2 cos (j))'^ + {VT 



yY + {^\ — x"^ cos 
^Jl — y"^ sin 0)2. Equating these two gives the consis- 
tency condition: 



xy 



so defining — \ — x"^ and / = 1 ~ y^ we have 



cos 



(1-/+)(1-/-) 



(B28) 



Now the Stokes vectors associated with v± (00) are 



5+(oo) 



^-(oo) 



1+ cos 29 
/+ sin 261 

A- \ 

cos 20 
/+sin20 
/ 



(B29) 



(B30) 



so the final Stokes vector of a state with initial Stokes 
vector 5*0 = S+{0) + aS-{0) is Sj = S+{oo) + aS-{oo): 



/+ cos 29 + al:^ cos 20 
/+ sin 29 + al:^ sin 20 

Thus the final polarization fraction, poo is: 



(B31) 



[1+ -al-f +AaI+I- cos'it 

[1+ + ai:;Y 

which after some simplification becomes: 
Poo = -F(X^cos20;;po) 



(B32) 



(B33) 



'1- 



4(1-P§)X^ 



2 ■ 



[(l + X2)-po(l-X2)cos2V'] 

where X ^ J7[0, 1) and -0 ^ U[0,2tt). In the simplest 
case where there is no initial polarization, po = 0, we 
have 



29 



which has mean value 



Poo = / dX = J - 1 « 0.57. (B34) 



/o 1 + ^^ 2 
More generally 



1 /■^'^ /"^ 
PaciPo) = TT / dXF(X2,cos2a;po)- 
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Pooipo) is a monotonically increasing function of po and 
increases from 7r/2 — 1 to 1 as po goes from to 1. 



APPENDIX C: ESTIMATING BL/2M AND 
CONFIDENCE INTERVALS 



In this appendix we provide details of how esti- 
mates and confidence intervals for the properties of any 
chameleon-like field can be extracted from measurements 
of the Stokes' parameters, I^, U and Q of a single object. 
We suppose that, in the absence of any chameleon field, 
the polarization angle of a given object is roughly in- 
dependent of wavelength in some interesting part of the 
spectrum (e.g. for UV to visible light). Since chameleonic 
effects die off as where A is the wavelength of light, 

we can roughly check this assumption by ensuring that 
the polarization angle is roughly wavelength independent 
for the larger wavelengths that are measured. We found 
in i jlVl above that the chameleon induced contributions to 
the expected Stokes' vectors oscillate fairly strongly with 
wavelength up until some critical oscillation wavelength 
Aosc- In all cases, we expect Aosc ^ 0{k). In addition 
to Aosc, there is another critical wavelength Acdt- Below 
Acritj the mean magnitude of the chameleonic polariza- 
tion signal is roughly independent of wavelength, whereas 
for A > Acrit ^ Aosc J the chameleon signal behaves as 
l/A^. 

We suppose that we have iVp measurements of the re- 
duced Stokes' parameter, Q/I-y and U/I-y, for a given 
object; we denote these measurements qt and ui respec- 
tively. We also require, for this analysis, that A > Aosc 
for all the measurements, and that any intrinsic (i.e. 
chameleonic) polarization be small i.e. ^ 100%. We 
define 5\ to be the spectral resolution of the measure- 
ments. In the weak mixing limit, to leading order, we 
have that the chameleonic contributions to Q/I-y and 
U / are given by: 



'Zcham 



• 2 A n 

-^o^EE'^«^(2A(n-r)) (CI) 
sin((5A(n — r)) 



■ cos(6l„ -I- Or), 



^cham 



5A(n - r) 

-^o^EE'^°^(2A(n-r)) (C2) 



sin((5A(n — r)) 

— TTT ^ — sm( 

i5A(n - r) ^ 



where Pq = {BL/2M) ; L is the coherence length of 
the magnetic field and B is its strength. N is the to- 
tal number of magnetic regions passed through, and M 
parametrises the strength of the chameleon to photon 



coupling. A = 7rA/2Acrit where Acrit = 47r 



"-eff 



m 



off 



rriu 



SA = 7r(5A/2 Acrit- The e„ define the 
angle of the magnetic field in the n*^ region relative 
to the direction of the light beam. Without any other 
prior information, we assume that these angles are es- 
sentially random. Now the total Stokes' parameters are 
Q — Qo + qcham and u = uq + Ucham- We assume that 
the non-chameleonic polarizations uq and qq depend on 
wavelength, but that, compared to the chameleonic con- 
tribution, they vary slowly. This will generally be the 
case, for instance, if both uq, go have a wavelength de- 
pendence similar to the the Serkowski polarization law 
[88j expected for polarization due to interstellar dust, 
i.e. uqPo oc exp(— ii'ln^(Amax/A)), for some K and Amax 
which we do not require to be the same for both uq and 
qq. Typically A 

max 

eOOOA and K w 1.15. We can then 
remove much of any intrinsic signal by simply smoothing 
the data over a scale on which uq and qo are expected 
to be fairly flat, to give and u^, and then subtracting 
this smoothed data from the original data. We define 
q = q — q^ and u — u — u^. 

We define yt for y^, with standard error <Ti, made at 
wavelengths Aj as follows: 

• We define some smoothing wavelength scale 
Asmooth and for each i define the Ji = 

{j ■ 2 |Ai — Aj)| < Asmooth}- 

• Ni is the number of elements in Jj. 

• We define = - S^^ T^jeJ, Vj/f^]- 

• Assuming that the yi are independent and dis- 
tributed N{pi, af) for some /j.^, we find the 5yi have 
standard error di ~ \/ — ^1 Si. 

We now have qi and from which the vast majority 
of any intrinsic signal should have been removed. We 
assume that any remaining intrinsic signal is sufficiently 
small compared with the noise as to be negligible. We 
check the accuracy of this smoothing process by simula- 



tions below. We define z\'^^ = qi/cfi and z-"^ = u^Jai. 

The chameleonic contributions to zf^ and are pre- 
dicted be and /3/i-"'' respectively where: 



,(<?) 



— E ^ikXk, 
fc=l 



(C3) 
(C4) 



fc=i 
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where (3 = Pq/2, h^k = H^k/cTi, Hik = Hu 



Htk and 



/^. = ^^^^^|^cos(2fcA.)5^, (C5) 



A 



with Ai = 7rAi/2Acrit and (5 A = 7r5A/2Acrit- We have also 
defined 



V2 



V2 

y'N -k 



cos(eW), 



r=0 
N-k-1 



(C6) 



(C7) 



Now if the Zi = yi/ci are just random noise with mean 
and variance A^ then u = u„ and: 



ijkl 



p'x'J2QkiMpL (CIO) 



= A^trM^^M^ - /) 
= -XhrV{l3), 

where V{P) — M^^ — I. Thus if there is only random 
noise we define E(^) — lq°^^°{(3) and we have: 



where el*"^ = Or+k + dr- When TV - fc > 1, and 
Yk are well approximated as independent identically dis- 
tributed N{0, 1) random variables. Since we assume that 
1 and the largest values of Hik occur for TV — A; » 1, 
we approximate the Xk and Yk as being independent and 
drawn from a A^(0, 1) distribution. The likelihood of find- 
ing Xk = Xk is therefore oc exp(— X|/2). 

Thus the probability density function with measure- 
ments Zj-*'' given /3 is (up to an overall Xk and f3 inde- 
pendent number Co): 



12k 



Defining the symmetric matrix Q thus Qik — X)i huhik, 
^1 + and Vk = J2i htkz['^^ we have: 

Since the first term is independent of both (3 and 
X we can incorporate it into a redefinition of the (3 
and Xk independent number Co i.e. Co Ci = 
Coexp(-^,4'^^V2). By defining Xk = ^p{X + 
PMj^^v), we have the new probability density, /, in terms 
of P, Zi and Xk'- 



fiP) = Do 



y/detMf3 ' 



where Dq is independent of (3. The /3 dependent term is 
now also independent of the Xk- We therefore define the 
likehhood of f3 given the qi data, i.e. the Vk, to be: 



L,(/3)=L,(0)- 



i/detM^ ' 



(C8) 



where L{Q) is the value of L when /? = 0. We define 
lg{f3)^ log L,{(3)/L,iO): 



2 



lq{l3) = '—v^M^^^v - 



■log detM;3. 



(C9) 



if'^P) = 2 l°g(^ + 7'(/3)) - XhxV{f3)) (CU) 



< -(l-A>7'(/3). 



(C12) 



with equality when (3 = and hence V{f3) = 0; generally 
ivV{l3) < with equality when /3 = 0. Thus if A = 1, 
which we should expect if the error estimates for the yi 
are accurate we have < for /? > 0. A more 

conservative approach would therefore be to use the data 
to check whether the scatter in the data points is as one 
would expect given the quoted errors, and if it is not 
extend the errors bars. We outline the method we use to 
do this in fJCT] below. Essentially, the highest frequency 
modes of any chameleonic signal, i.e. those with k N — 
1, also produce the smallest contribution to the overall 
signal; all other modes are approximately constant over 
wavelength scales of about Acrit/(-^— !)• Thus, provided 
there are enough data points, we can use the variance 
of the data points on scales < Acrit/(^ — 1) to estimate 
their error. 

We make a similar set of definitions for the Ui data, 
for which Z„(/3) is the log-likelihood, and define the total 
log-fikelihood to be ?(/3) = lu{(3) + lq{f3). We define the 
maximum likelihood estimate of [3, [3, to be the value 
of 13 which maximises 1{I3). There are now a number 
of approaches we may take to estimate the confidence 
intervals. The simplest approach is to assume that the 
values of (3 are normally distributed with some variance 



CT^. We then estimate tr^ as: 



4 



IMP) 



(C13) 



A 95% confidence interval for (3 is then estimated to 
he (3 = (3 ± 1.96(7^. Since the quantity we are ac- 
tually interested in is a; = \B\L/2M = \fWi we dis- 
play all confidence limits as constraints on the value of 
X. We refer to this as the normal approximation and 
label it (NA). The second approach is to assume that 
r{(3,(3) = 2(10) - l{(3)) - xl, which should hold as the 
number of observations tends to infinity; the approxi- 
mate 95% confidence interval for /3 is all /3 for which 
r{(3, f3) < 3.84. We transform this into an approximate 
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confidence interval for x by taking x — y 2/3. We refer 
to this as the approximation, labeUed (x^)- A more 
robust approach to estimating confidence intervals is to 
bootstrap the Ui and qi data. We make B bootstrap data 
sets constructed by resampling with replacement Np data 
points for both Ui and the cji from the original data. For 
each bootstrap data set we construct the /3 and in 
the same way as was done in the normal approximation, 
defining them to be: /3* and ^ . Now there are a number 
of different bootstrap methods for estimating confidence 
intervals. We use the bootstrap-t method which generally 
has better convergence than the usual bootstrap method. 
We assume that the distribution of i = (/3 — /3) /ct/j is well 
approximated by the distribution of the the bootstrap 
parameter t^, — {(3^ — j3)/ap^,. The lower limit, (3^ of 
the bootstrap-t 100(1 — 2a)% confidence interval for j3 is 
therefore given by: 

B-B P -P 
a = P{t^J± < '=^J1) ^ Pit > -t) = l-P{t < -tj. 

where = We estimate P{t < —t^) by P*(t* < 

-to) = G Boot {-ta) where 



confidence limits for x = BL/2M: 



<s) = GBoot(/3) = 



< s) 
B 



Thus we have: 

Similarly the upper limit is: 

We define the central estimate of /3 to be /3„ 



^/3^Boot( 1/2)- If P is an unbiased estimator for P then 
GgQQj(l/2) = and we have Pm = P. We label this ap- 
proximation (Bt). In all cases given below we have used 
B = 5x 10** bootstrap resamplings. When the error bars 
are rescaled as described above and in 1\ we construct 
confidence intervals in both the normal approximation, 
and using the technique; we label these approaches 
(NA - a) and (x — o") respectively. 

Using starlight polarization data for three objects from 
the WUPPE spectrograph with a spectral resolution of 
leA and a spacing between data points of 2A, we are 
able to find useful constraints. In all cases we take 
Acrit = 608A. A more thorough analysis would attempt 
to also fix Acrit, given some reasonable priors about the 
electron density, Ue, the chameleon mass, m^, and the co- 
herence length. In all cases we take Asmooth = lOOA. We 
find that for 75A< Aginooth ^ 200A our results do not 
depend greatly on Agmooth- For the first star, HD2905 
{d = 880 pc), we have the following approximate 95% 



(6.36;°;^2) X 10-2 
(6.36;j;°?) X 10-2 
(4.68;i;^4) X 10-2 
X - (4.1lt°:?^) X 10-2 



X 
X 
X 



(NA), 

ix% 

(Bt), 

(NA-a), 
(x'-a). (C14) 



Even when the error bars are extended as described 
above, we have that r{x) = 2l{x) = 91.9; indicating that 
the maximum likelihood estimate for x deviates from 
by more than 9.5a in the ~ ^ approximation. We note 
that for HD2905, Acrit ~ 610Ais a local maximum of the 
likelihood the MLE estimate for x. 

The last three techniques are all in rough agreement. 
We see the same behaviour in the analysis of simulated 
data that we have undertaken; these simulations also 
show that the last three techniques are most accurate, 
and are robust to the actual errors being larger than the 
quoted ones. From these simulations we also find that 
by approximating the Xk and Yk as independent identi- 
cally distributed N{0, 1) random variables, we reduce the 
likelihood of the MLE for x, but do not greatly alter the 
value of the MLE. 

For two other objects (again assuming Acrit — 608A), 
we find similar results. For IID39703, at d = 880 pc, we 
find: 



X 



X = (8.69+'-''^'^ - i"-2 



X — 



-1.39) X 
.63^ 



{7.59tH^) X 10-2 



x=(8.iii}:^,°)xio-2 

X = (8.1lli; 



-ID X 10-2 



(NA), 
(Bt), 

(NA-a), 



(C15) 



When the error bars are extended as described above 
we have r{x) — 2l{x) — 74.9. This implies that, in the 
x2 — (7 approximation, a; = is more than 8.6f7 from 
the maximum likelihood estimate of x. For HD34078 
{d — 610 pc) we have: 



= (9.95t^:I?) 


X 10-2 


(NA), 


- (9.951?:°^: 


) X 10-2 


ix'), 


= (8.58i?-i^: 


) X 10-2 


(Bt), 


= (9.4it^:^^) 


X 10-2 


(NA-a), 


- (9.4iiig^: 


) X 10-2 


ix'^a). 



(C16) 



In this case, r{x) — 2l{x) — 84.9 when the error bars are 
extended. This again implies that the maximum likeli- 
hood estimate for x deviates from by more than 9a in 
the x^ approximation. As expected from simulations, 
the last three techniques are all in rough agreement. 

If we assume that the same value of BL/2M should 
be appropriate for all three objects (which may not nec- 
essarily be the case), then combining all three data sets 
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we find the following 95% confidence intervals: 



1. Extending the Estimated Errors 



X = (7.95tgJi) X 10-^ 
X = (7-95t°;i) X 10- 
X = (6.25t J:23) X 10- 



(6.03+° 



+0.76\ 
87/ 



X 10" 



X = (6.03tg;^^) X 10- 



(NA), 
(Bt), 

(NA-a), 
ix'-a). (C17) 



The log-likelihood of the MLE of x when the errors have 
been rescaled is r = 21ogZ(i;) = 214; indicating a more 
than 14.6cr deviation from in the approximation. 
As with all three objects separately, wc sec that there is 
rough agreement between the last three approaches, al- 
though, as was the case for the three objects separately, 
the error bars are widest in the bootstrap-t approxima- 
tion. 

Combining all the data in the standard approach (as- 
suming the same value of a; = BL/2M is appropriate 
for all), and using the bootstrap-t method, we find the 
following 99.9% confidence intervals: 



^ = (6.25-%) X 10-, 



(99.9%). (C18) 



Using the bootstrap-t method, inasmuch as the resolu- 
tion of the bootstrap distribution allows, we find that, 
defining ax* — ap^/x*, the distribution of 6*0(2;* — 
x)/ax* + Si is approximately A^(0, 1) for some 5*0 and 
Si . If we assume that the distribution of (i* — x)/ax* is 
a good approximation to that of {x — x)/ax, we have: 

^ = ^ = (6-27 ± 0.58) X 10-2, (C19) 

where this time the quoted error bars are Icr. This cor- 
responds to a more than lO.Tcr deviation from 0, and 
provides the following 95% and 99.9% approximate con- 
fidence intervals: 



— = (6.27 ± 1.14) X 10-2, 
2M ^ ' ' 

BL 



(95%), (C20) 



^ = (6.27 ± 1.91) X 10-2, (99.9%). (C21) 

Using our estimated values for B and L we have at 
99.9% confidence: 

M = (1.471^:^^) X lO^GeV, (99.9%). (C22) 

Although this analysis is only preliminary, it does 
appear as if there is a reasonably significant, and ro- 
bust, statistical preference towards the existence of a 
chamclcon-likc field in the starlight polarization data of 
the three objects we have considered here. A fuller analy- 
sis would have to take into account more, even all, compa- 
rable starlight polarization measurements. Additionally 
one would also wish to fit for Acrit- 



In this subsection we provide further details of how we 

extend the errors bars on the data to better mask the 
observed small scale scatter. We expect any chameleon 
induced fluctuations of the polarization on wavelength 
scales smaller than Acrit /(A^ — 1) to be small compared 
with that on larger scales between Acrit/(A-- 1) and Acrit. 
In all cases we estimate Acrit/(A^ ~ 1) ^ 16A. For each 
smoothed data point {Xi,Ui,qi), we use the data points 
labelled j with 2|Aj — Aj | < 16A to estimate the random 
(or non-chamclconic) scatter in the data. The estimated 
standard errors in the smoothed data points are tTj. We 
define, as we did above, Ji = {j ■ 2 |Ai — A^ l < Asmooth}, 
where this time Agrnootii = 16 A. 

We assume that the data points in Ji have mean /i and 

standard error yoJ+cJcj^, where da'^ is to be found (/U 

and a will be different for the qi and the tt,). For data 

points Xj , with estimated standard error aj , where j G Ji , 
we estimate /z by /x its maximum likelihood estimator: 



(C23) 



Similarly for each i, we estimate Sa'^ by its MLE daf 
which satisfies: 



E 



= E 



2 ^a^+5a^- 



(C24) 



If no solutions to this equation exist, then wc take Saf = 
0. Finally wc smooth the Saf over a IGOAsmoothing 
scale, giving 6af. We take the final enhanced error to be 

= \/ o'i + Saf. Although this procedure is rather ad 
hoc, by enhancing the error bars, we err on the side of 
caution and reduce the probability that under-estimated 
error bars result in a spurious detection of /3 ^ 0. 

We present estimated confidence intervals where the 
error bars have been extended using the normal approx- 
imation and the approximation; we label these two 
approaches (NA-cr) and (x^ — a) respectively. 



2. Estimating upper bounds on BL/2M 

It is also possible to find upper confidence limits on 
BL/2M, simply from the observation that the com- 
ponent of polarization perpendicular to the mean po- 
larization angle is smaller than some upper bound i.e. 
i-PL(A)| < Pmax- Suppose the observations of P± are Qi, 
and that we have the maximum value of the qf < p^ax- 
If a chameleon field is present, and assuming that the 
polarization angle of any intrinsic polarization is roughly 
constant, we predict = fPi^k-^khkif ■ If the intrin- 
sic polarization angle is not constant then we will gen- 
erally be biased in favour of larger values (3, and so this 
approach can also be trusted to provide upper bounds 
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on BL/2M. Using numerical simulations we can esti- 
mate the distribution of w = max(^^ X^j/i^i)^ and for 
< a < 1 calculate Wa the probability that: 

P{w < w^) = a = P{f3^w < (3^Wa). 

We defining f/j^^^x — uiSixqf we then have: 



Thus qiia^/wa = Pmax/'"'" = is a estimate of the 
100(1 — a)% upper confidence limit on Generally this 
is an over-estimate of the true upper confidence limit, and 
so (3 < f3a with at least 100(1 — a)% confidence. 



< 9max/w^«) = 1 - a. 
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